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The ionization of single cosmic-ray particles capable of penetrating 12 cm of lead has been 
measured with an ionization chamber and a linear amplifier. Measurements made at 14.7 
atmospheres with argon indicate an average specific ionization of 71 ion pairs per cm at N.T.P. 
The results apply on the whole to mesotrons with an energy greater than 2X 10* ev. The specific 
ionization of maximum frequency of occurrence is approximately 67 ion pairs per cm in argon. 
The average value for the ionization obtained with the ionization chamber at 14.7 atmospheres 
appears to be about 10 percent less than the total ionization. The average value obtained 
corresponds to the value 48 ion pairs per cm in air at N.T.P. as measured by a method that 
includes only collisions involving 10* ev or less. This result is in good agreement with the cloud- 
chamber results of Corson and Brode and of Hazen. 


I. INTRODUCTION 


LTHOUGH the ionization chamber has been 
used very extensively for the measurement 
of ionization currents and for the measurement of 
ionization pulses of large numbers of cosmic-ray 
particles such as form a shower or burst, the only 
investigator previous to the present work who has 
used the instrument for the measurement of the 
ionization of single cosmic-ray particles is W. F. 
G. Swann.' Swann used a vertical cylinder filled 
with argon to a few atmospheres pressure. No 
counter control was used, and all types of pulse 
were recorded. A rather elaborate analysis had to 
| * This is a condensation ( ther with added 
material) re thesis presented in yey 1942 to the 
Graduate Division of the University of California in 
_— fulfillment of the requirements for the degree of 
** Now at the Western Regional Research Laboratory, 
Albany, California. 
‘W. F. G. Swann, Phys. Rev. 43, 961 (1933). 
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be made to obtain the probable distribution of 
path lengths inside the cylinder. In abstract form 
Swann?’ has also mentioned use of counter control 
with the ionization chamber. 

The work reported in the present paper was 
begun with the purpose of obtaining a new and 
independent value for the specific ionization of 
fast mesotrons, and if possible, of developing a 
method for obtaining accurately the specific 
ionization of a single mesotron. The momentum 
could be measured at the same time by means of 
a cloud-chamber measurement of the curvature 
of the track in a magnetic field. Although both 
the ionization and momentum can be obtained 
from cloud-chamber photographs, the diffuse 
tracks needed for accurate droplet counts make 
determination of the momentum inaccurate if 
only a single cloud chamber is used. 


2 W. F. G. Swann, Phys. Rev. 45, 258 (1934). 
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Il. APPARATUS 
A. General Description 


A schematic diagram of the apparatus is shown 
in Fig. 1. The pulse produced by the motion of 


YEO 


Fic. 1. Schematic diagram of ap tus. The chamber, 
counters, and 38 tube are mounted inside a shielded box, 
which is hung by a single piano wire from rubber cushioned 
supports to minimize microphonic effects. 


the ions liberated by the primary mesotron or 
other particle was detected at the grid of the 
38 tube. A high ‘electrostatic field was used to 
collect the ions. The pulse was pre-amplified by 
the 38 tube and further amplified by the linear 
amplifier, which delivered the amplified pulse to 
the vertical deflection plates of the oscilloscope. 
The electron beam of the cathode-ray tube was 
deflected off the screen until a coincidence pulse 
caused the sweep-circuit’ thyratron to break 
down. The thyratron breakdown produced a 
single trace upon the screen. If the coincidence 
pulse was accompanied by an ionization pulse, as 
was ordinarily the case, the pulse appeared upon 
the screen at the start of the trace. Because of the 
difficulty of obtaining accurate visual observa- 
tions of the length of the pulses, and because of 
the low counting rate—one in ten to twelve 
minutes—the apparatus was made self-recording 
and nearly automatic. In the arrangement finally 
adopted, the traces were photographed by means 
of a camera equipped with an automatic rewind, 
actuated by the pulse from the coincidence 
counter set. A cam system was used to cause the 
camera gear to be turned the proper fraction of a 
revolution after each pulse. The apparatus ran 
automatically for eight or nine hours on one 
loading of film. 


B. The Ionization Chamber 


A schematic diagram of the chamber in cross 
section is shown in Fig. 2. The chamber was 
cylindrical. The walls were designed to withstand 
a pressure of at least 70 atmospheres, although 
for the experiments reported here much lower 
pressures were used. Three plates were used in 
order to secure flexibility of operation. Thus, the 
two outer plates could be connected independ- 
ently to high voltage sources of the same or 
opposite sign. In some cases it is desirable to use 
voltages of opposite sign so that for particles that 
produce equal ionizations in the two halves of the 
chamber no net effect is observed at the center 
electrode. The center collecting plate was insu- 
lated by use of amber that was coated with 
ceresin wax to reduce surface conductivity. The 
spark plugs used as leads for the high voltage 
plates were plugged with wax to make them gas- 
tight. Hard rubber insulators were found to be 
satisfactory supports for the high voltage plates 
at voltages up to 7500. 

Microphonics arising especially from the vibra- 
tions of the plates were a serious problem 
throughout the course of the work. All three 
plates were made 3%” thick to minimize these 
effects, and the whole apparatus was hung by a 
single piano wire from rubber-cushioned sup- 
ports. These measures reduced the sensitivity of 
the system to vibration to the extent that only a 
very few pulse records had to be discarded be- 
cause of the presence of microphonic disturbances. 


C. The Electrical Circuits 


The four-stage linear amplifier was of the 
negative-feedback type; the feedback was pro- 
vided by use of un-bypassed cathode resistors 
and by two feed-back loops. The circuit was built 
up according to the design of Waddel® and 


differed from that of Waddel in only minor de- - 


tails. The amplification of the amplifier itself, 
including the 38 tube but not the oscilloscope 
amplifiers, was approximately 60,000. The cali- 


bration of the amplification system is discussed . 


below. 
The double-coincidence counter set was of 


standard design, and will not be described com- 
pletely. The counter tubes were quenched by 


*R. Waddel, Rev. Sci. Inst. 10, 311 (1939). 
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means of the Neher-Harper‘ arrangement, and 
the coincidences were detected by use of the 
Rossi® coincidence circuit. The resolving time of 
the counter set was about 2X 10~ second. A time 
delay section employing an 885 thyratron was 
used to delay closing of the rewind-relay contacts 
until after the ionization pulse had been photo- 
graphed. This measure was required because of 
the relative slowness of the ionization pulse com- 
pared to the coincidence pulse and because 
closing of the relay contacts during recording of 
the pulse would have caused sufficient sparking 


" to obliterate the trace of the pulse. In order to 


operate the camera motor from the very short 
pulse from the coincidence counter, a standard 
self-holding telephone type relay with two 
energizing coils was used. The relay was closed by 
the coincidence pulse and remained closed until 
the holding current was stopped by the breaking 
of the relay contacts by the cam on the motor 
shaft. 

The high electrostatic field required for the 
collection of the ions was obtained by use of a 
voltage doubler circuit that provided either plus 


or minus 2900 volts with respect to ground, or 


both. When it appeared desirable to use even 
higher voltages, a neon sign transformer was used 
in the rectifier set, and the circuit was changed so 
as to provide a single voltage variable from 0 to 
7500 volts below ground. No results obtained by 
use of this voltage, however, are presented here. 
Since only small currents were drawn from the 
power supply of the high voltage unit, the little 
filtering required was provided by two sets of 
high resistances and one-microfarad condensers. 
There remained, however, some 60-cycle voltage 
in the background. This ripple apparently arose 
from pick-up of atmospheric fields, and was only 
partially removed by careful grounding of the 
shielding conductors. 


Ill. CALIBRATION OF THE AMPLIFYING SYSTEM 


The over-all amplitude and frequency charac- 
teristics of the amplifying system, including the 
oscilloscope amplifiers, were measured by use of 
alternating voltages obtained from a signal gen- 
erator. From the frequency characteristic it was 
established that the amplification of the system 

*V. H. Neher and W. W. Harper, Phys. Rev. 49, 940 


(1936). 
5B. Rossi, Nature 125, 636 (1930). 


was constant for frequencies corresponding to all 
collecting times used in the experiments.* From 
the amplitude characteristic it was verified that 
the length of the pulses was proportional to the 
voltage produced upon the grid of the 38 tube. 

The amplitude and frequency characteristics 
mentioned above were only relative. The system 
was calibrated for what might be called the 
“absolute sensitivity” by use of a polonium 
alpha-particle source within the ionization cham- 
ber. The chief requirements for a good calibration 
are: (1) The pulses obtained must be of the same 
order of magnitude as the cosmic-ray pulses 
studied. (2) If possible, the particles should trav- 
erse the chamber in the same way as the cosmic- 
ray particles are expected to do. (3) The probable 
error of the result must be made small by obser- 
vation of a large number of alpha-particle pulses. 
The method of calibration finally adopted in- 
volved use of an alpha-particle gun that confined 
the beam within a narrow solid angle in a direc- 
tion parallel to the plates, and close to the voltage 
plate. The chamber was at a reduced pressure 
(~10 mm Hg) and the range of the particles was 
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Fic. 2. Schematic cross-section diagram of the 
ionization chamber. 


great enough to enable the alpha-particles to 
traverse the length of the plates much as a 


* Since these results were obtained with low impedance 
sources, they do not include an effect to be discussed below, 
namely, a variation of sensitivity with collecting time 
- gd the leakage of charge from the grid circuit of the 

tube. 
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cosmic-ray particle would have done. By use of 
the ionization versus range data given by 
Rutherford, Chadwick, and Ellis,’ the total num- 
ber of ions expected in the sensitive region was 
calculated by graphical integration of the ioniza- 
tion. The resulting value of the sensitivity of the 
amplifying system was 22,300 ions per inch on 
the screen at the oscillograph gain setting used 
for the majority of the measurements. The 
standard deviation for this result is 2.0 percent. 

One of the chief difficulties of the experiments 
was the fact that the time of collection of the 
ions, for high pressure, was not small compared to 
the time constant of the grid circuit of the 38 
tube. This time constant is effectively the product 
of the leakage resistance from grid to ground and 
the capacity of the grid-collector system with 
respect to ground. The time constant for the 
apparatus used here was about 0.01 sec. The 
variation of pulse size with collecting time (or 
collecting field) was determined by use of the 
alpha-particle method mentioned above. By com- 
parison of the resulting curve with the saturation 
curve for alpha-particles given in Jaffe’s paper on 
columnar ionization,’ it was found that in order 
to have the uncertainty caused by the collecting 
time effect less than 10 percent, the collecting 
time must be less than 0.005 sec. The collecting 
time must be less than 0.003 sec. for an uncer- 
tainty less than 5 percent. In addition, each value 
obtained at 27 atmospheres was increased by 25 
percent, at 21.5 atmospheres by 12 percent, and 
at 14.7 atmospheres by 5 percent, to correct for 
loss of pulse size from the collecting time effect. 
The values for the collecting time used above 
were determined for the present calculations by 
use of the value 1.7 cm?/volt sec. for the mobility 
of argon containing slight impurities of oxygen.° 
(Argon of 96 percent purity was used in this 
work.) 


IV. SOME FURTHER CONSIDERATIONS 
AND ASSUMPTIONS 


In the calculation of the specific ionization of 
the penetrating particles from the observed pulse 


7 Rutherford, Chadwick, and Ellis, Radiations from 
Radioactive Substances (Cambridge University Press, 1930), 


. 80. 
4 8G. Jaffe, Ann. d. Physik 47, 303 (1913). 
*L. B. Loeb, Kinetic of Gases (McGraw-Hill Book 
Company, Inc., New York, 1934), p. 619. 
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sizes, the following assumptions were made: 
(1) Recombination can be neglected. Clay!® has 
shown that for air in a shielded ionization cham- 
ber, with fields of the order of 2000 volts/cm and 
pressures of 15 to 25 atmospheres, the lack of 
saturation is only one or two percent. Since the 
pressures used in this work varied from 14.7 to 27 
atmospheres, and the collecting field for most of 
the measurements was 1900 volts/cm, recombi- 
nation can reasonably be neglected as a source 
of error. (2) The measured ionization is pro- 
portional to the pressure. Clay" has shown that 
the measured ionization in argon is proportional 
to the pressure, to within a few percent up to 60 
atmospheres. 

The path length used in the calculation of the 
probable values of the specific ionization is 17.0 
cm. This length is not exactly that of the plates, 
but is a corrected value slightly greater than the 
length, 15.9 cm. The corrections take into ac- 
count the effect of the fringing fields at the ends 
of the plates and also the fact that the Geiger 
counters subtended an appreciable solid angle. 
The counters were disposed so as to cover the 
entire region between two of the plates. 

By use of the alpha-particle calibration and the 
above value of the probable path length, the 
“probable value of the specific ionization’’ was 
calculated, and it is equal to the observed ioniza- 
tion (corrected for collecting time effect) divided 
by the path length and by the pressure in atmos- 
pheres, reduced to standard conditions. 

The expression ‘‘probable value of the specific 
ionization”’ should be used to denote the results, 
since it is not certain that a given pulse was pro- 
duced by a single particle nor that the particle, 
assuming it to be single, traversed the chamber in 
the expected way. On the whole, however, the 
results can be ascribed to the action of single 
penetrating particles, presumably mesotrons. 
Starr” and others" have shown that the great 
majority, of the order of 90 percent, of double 
coincidences obtained with a few centimeters of 
lead between the counters is caused by single 
penetrating particles that do not produce showers. 
In the present experiments, the walls of the 


10]. Clay and K. Oosthuizen, Physica 4, 527 (1937). 
11]. Clay and M. Kweiser, Physica 5, 725 (1938). 
12M. Starr, Phys. Rev. 53, 6 (1938). 

13S. H. Neddermeyer and C. D. Anderson, Phys. Rev. 


51, 884 (1937). 
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chamber, the bloek of lead, and the concrete roof 
above the apparatus are equivalent as far as 
jonization loss is concerned, to about 12 cm of 
lead. Using the Bethe-Bloch value of 1.5X10° 
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Fic. 3. Pulse-size-frequency distributions, plotted as 
block diagrams. The abscissa is the “probable _——— 
jionization,”’ that is, the specific ionization to be ascribed to 
the pulse assuming it to be produced by a single particle 
that traverses the chamber with path length 17.0 cm. 


ev cm*/g, for the ionization loss of a fast mesotron 
at the minimum of the ionization versus energy 
curve, we see that the results apply on the whole 
to mesotrons with energies above 2X 108 ev. 


V. RESULTS AND DISCUSSION 


Measurements were made at pressures of 27.0, 
21.5, and 14.7 atmospheres of argon. Higher 
pressures were not used because of the limitations 
on the collecting time mentioned in Section III. 
About 1000 pulses were obtained in the form of 
' photographic records, of which 700 were suitable 
for measurement. The block diagrams of Fig. 3 
show the frequency of occurrence of certain 
specific ionization intervals of widths indicated. 
From the study previously mentioned of the 
variation of sensitivity of the amplifying system 
with time of collection, it was found that at 22 
atmospheres there was an uncertainty of 20 
percent in each individual result, 10 percent at 
21.5 atmospheres, and about 5 percent at 14.7 


atmospheres. In spite of this uncertainty, and of 
the previously mentioned corrections to the 
ionization values at the three pressures used, the 
maxima of the three distributions all appear to be 
in the same region. There is of course a pro- 
nounced broadening of the distributions at the 
higher pressures, which can be attributed to the 
collecting time effect. 

The frequency distribution curve is consider- 
ably broader than would be expected if single 
particles having an energy spectrum such as that 
found by Blackett and obeying the Bethe-Bloch 
ionization law were passing through the chamber 
in the directions defined by the two Geiger- 
Miller counters." In particular, it is to be noticed 
that there is a considerable number of pulses 
smaller than that most frequently occurring, as 
well as a large group of pulses longer than that 
most frequently occurring. 

This broadening cannot be attributed to the 
“collecting-time effect,’’ since the error in each 
result introduced by this effect at 14.7 atmos- 
pheres is only about 5 percent. There are, how- 
ever, additional factors peculiar to the experi- 
mental method that may be responsible for the 
observed breadth of the curve. Some of these are: 
(1) Single penetrating particles may pass through 
the counters without traversing the full length of 
the sensitive region of the chamber to produce a 
small pulse. The coincidence pulse would most 
likely be produced in this case by scattering of the 
particle in the walls of the chamber or by pro- 
duction of a knock-on electron at the outer edge 
of the chamber. (2) Showers may trip the counters 
and one or more of the particles of the shower 
may pass obliquely through the sensitive region 
of the chamber to produce a small pulse. (3) Acci- 
dental coincidences may be caused by two inde- 
pendent cosmic-ray particles, one of which by 
chance may also pass obliquely through the 
sensitive region of the chamber to produce a 
small pulse. It is expected that showers are mainly 
responsible for the existence of both the very 
small pulses and the very large pulses observed. 

The value taken for the specific ionization of 
maximum occurrence is 67 ion pairs per cm in 
argon at N.T.P. The average specific ionization is 


4 For information on the question of the constancy of 
specific ionization of cosmic ray particles with energy, see 
. E. Hazen, Phys. Rev. 65, 259 (1944). 
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_ also of interest. Using only the 14.7-atmosphere 
results we find by counting squares that the 
average specific ionization is 71 ion pairs per cm. 


In order to obtain maximum accuracy in - 


measurement of the most frequently occurring 
pulses, pulses with lengths more than about three 
times that most frequently occurring extended 
beyond the limits of the oscilloscope screen and 
were not measurable. By use of minor modifica- 
tions of the photographic technique, the range of 
investigation can be modified so as to include any 
range of ionization, from that of single particles 
to that of importance in the study of showers and 
bursts. 

If one uses Hopfield’s assumption’® that the 
ionization is proportional to the number of extra- 
nuclear electrons, then the specific ionization for 
air should be about 14.5/18 that for argon. For 
the average specific ionization in air from the 
results at 14.7 atmospheres, one obtains a value 
- of approximately 60 ion pairs per cm. 

Although a minimum pressure of 14.7 atmos- 
pheres was used in the chamber, it cannot be 


assumed that the average, ionization obtained in 


this work represents the total ionization.’® A 
secondary produced in the gas can impart only a 
limited energy to ionization and excitation within 
the sensitive region of the chamber. Assuming a 
constant ionization of 60 ion pairs per cm, and an 
average energy of 25 ev in argon for production 
of one ion pair,!” we find the maximum energy of 
dissipation of a secondary to be 3.75 X 105 ev, if 
we assume that the secondary passes through the 
sensitive region in a straight line with path length 
17.0 cm. This assumption is doubtful, however, 
since scattering of the particle would have the 
effect on the whole of shortening the path length 
by deflecting the particle into one of the plates. 
Secondaries with energies above 375 kev will be 
able to contribute only this same energy. Thus 
although it is rather difficult to calculate a pre- 
cise upper limit to the secondary energies that are 
detected, for a rough calculation we take 375 kev 
as an effective upper limit and consider that 


secondaries with energies up to 375 kev are com- . 


J. J. Hopfield, Phys. Rev. 43, 675 (1933). 
am indebted to Dr. Wayne E. Hazen for pointir this 
én - me and for valuable assistance in connection with the 


point 
17 See reference 7, p. $I. 
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pletely detected, and that secondaries with 
energies above 375 kev are completely undetected. 

The total collision loss @, in air is given by 
Rossi and Greisen'* from the Bethe-Bloch theory 
to be 1.8 X 10° ev cm?/g in air. At one atmosphere 
the total average energy loss is 2300.ev/cm. Rossi 
and Greisen also show that the average energy 
loss when only collisions up to 10‘ ev are con- 
sidered is 28.5 percent less than the total, or 
€10¢=1700 ev/cm. We now wish to add to the 
latter figure the energy loss caused by collisions 
to secondaries between 10‘ and 3.75 X 105 ev. The 
probability of production of a secondary electron 
of energy & by a fast mesotron of high energy is 
given fairly accurately for air by the Rutherford 
formula!® 

98dE ev 


x(E)dE= 
E2 
The average energy loss éz; caused by all such 
secondary electrons in the range FE, to Ez is 
Es E2 
f Ex(E)dE = 98 —. 
E 
For E,=10* ev, E2.=3.75X105 ev, 
98 In 37.552360 ev/cm. 


Since the ionization chamber effectively meas- 


ures we find a fraction 2060/2300 
=90 percent for the fraction of the total col- 
lision loss measured by the ionization chamber as 
used in the present work. Since the cloud- 
chamber result (for energies less than 10* ev) is 
28 percent below the total ionization, the cloud- 
chamber result is about 20 percent less than the 
ionization chamber result. Thus for comparison 
of the present ionization chamber results with 
the cloud-chamber results, we should use thé 
value 60— (0.260) =48 ions pairs per cm path 
in air at N.T.P. 

Corson and Brode*®® found a minimum ioniza- 
tion of fast electrons in air of about 50 ion pairs 
per cm. Bagley” found an average ionization of 


18B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 240 


(1941). 
19B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 243 
(1941). See also reference 14. 
= Ly “— and R. B. Brode, Phys. Rev. 53, 773 (1938). 
G. Bagley, Ph.D. Thesis, University of California 


(1941). 
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cosmic-ray particles in nitrogen of 43 ion pairs- 
per cm. Hazen” has recently obtained a value of 
50 ion pairs per cm for the average specific 
jonization in air of a large number of tracks of 
penetrating rays obtained at sea level and 100 
feet underground. 

It is of interest to compare the present results 
with those obtained by other workers using 
ionization chamber methods. Swann* obtained a 
value of 50 ion pairs per cm in argon by use of the 
uncontrolled ionization chamber having a thin 
wall and with little shielding around the chamber. 
Stuhlinger™ in work with a proportional counter 
at one-half atmosphere, found a value for the 
specific ionization of the hard component in air 
of 35 ion pairs per cm.*® Stuhlinger in this same 
work, and Swann?* in work reported only in 
abstract form, have found a second maximum in 
the ionization curve. Indeed Stuhlinger gives a 
curve showing six or eight fairly well-defined 
maxima which he ascribes to various combina- 
tions of primary hard particles and softer second- 
ary particles. The primaries are assumed from 
the position of the first maximum to have a 
specific ionization of about 35, and the second- 
aries then are found to have an ionization of 
about 50 ion pairs per cm, referred to air. 

There is some evidence of such a second maxi- 
mum in the 21.5-atmosphere results, but none in 
either the 27-atmosphere nor in the 14.7-atmos- 
phere results. Since the low pressure results are 


the most reliable, it must be concluded that the. 


present results indicate that there is only one 
maximum in the frequency distribution. 


VI. CONCLUSION 


The results obtained in this paper for the 
average value and for the most frequently oc- 


= W. E. Hazen, Phys. Rev. 65, 259 (1944). 

%W. F. G. Swann, Phys. Rev. 43, 961 (1933). See also 
footnote, T. H. Johnson, Rev. Mod. Phys. 10, 211 (1938). 

*E. Stuhlinger, Zeits. f. Physik 108, 444 (1942). 

* Information has just been received of the publication 
of a paper “On the use of the linear amplifier for the 
measurement of the ionization of single particles’’ by S. A. 
Wytzes and G. J. van der Maas, Physica 10, 419 (1943). 


No details, however, are available at the present writing. 
_ %*W.F. G. Swann, Phys. Rev. 61, 393 (i 


942). 
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curring value of the specific ionization obtained 
under the described conditions are considered to 
be accurate to within eight or ten percent. The 
principal sources of error in these quantities are 
(1) the error of calibration of the sensitivity of 
the amplifying system, and (2) the fact that the 
pulses were classified into only a relatively few 
large groups. This classification was dependent 
upon the diffuseness of the trace, which in turn 
was determined by the background disturbances 
in the 38 tube. 

In evaluating the question of the accuracy 
with which a single particle can be measured, we 
are probably upon safer ground if we examine the 
breadth of the experimental curve of Fig. 3 than 
if we rely upon estimates of the uncertainties that 
may be supposed to be introduced into each 
measurement. ‘The curve has a _ half-breadth 
which appears to be some 25 or 30 percent of the 
most frequently occurring value. Although there 
is no reason to suppose that such a large error as 
25 percent was inherent in the method of measure- 
ment of the ionization in the chamber, neverthe- 
less, the presence of knock-on electrons from the 
walls or of pulses due to showers leads to essen- 
tially the same effect. Although use of the cham- 
ber in conjunction with a cloud chamber might 
remove the uncertainty as to whether a given 
pulse is associated with an air shower, there 
appears to be no way of being sure that a knock- 
on electron is not produced in the wall of the 
chamber to give a misleading impression of the 
specific ionization of the particle. Thus the 
usefulness of the method outlined in the present 
paper is greatest for work in which the most 
frequently occurring value, based upon as large a 
number of pulses as possible, can be used as the 
measure of the ionization. 
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A Study of Time Variations in the Cosmic-Ray Directional Intensity Distribution* 
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Evidence of time variations in the cosmic-ray directional intensity distribution has been 
obtained, by a statistical analysis of directional intensity measurements carried out at St. Louis 
in 1943, The variations, occurring within periods of twelve to twenty-eight hours, seem most 
pronounced near zenith angles of +20° and +70° in the east-west plane, and are not found in 
the north-south plane. These observations would appear to be explained by the effect of 
terrestrial magnetic fluctuations in modifying the allowed cones of certain primary cosmic-ray 
particles, these particles having energies grouped near 0.2 Stérmer and 0.6 Stérmer. If this 
hypothesis is correct, the east-west positional symmetry of these variations constitutes new evi- 
dence, supporting previous evidence of Schremp, that positive and negative particles with the 
same ¢/m exist in the primary cosmic radiation. Earlier results of Ribner, Cooper, and Schremp, 
on the static fine structure pattern of cosmic-ray directional intensity, are also supported by the 
present data. The time variations in the intensity distribution reported here are variations in 
the fine structure pattern, but appear to be independent of the static prominences previously 
reported. The method for rapid exploration of the cosmic-ray directional intensity distribution 
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is described. The multidirectional cosmic-ray telescope used is discussed in some detail. 


I. INTRODUCTION 


ROM theoretical considerations, based on 
the theory of the allowed cone! and the 
theory of atmospheric absorption of cosmic rays, 
Schremp predicted*? that a fine structure pattern 
should exist in the cosmic-ray directional in- 
tensity distribution, and indicated* how a study 
of this pattern should provide information con- 
cerning the nature, energy spectra, and absorptive 


properties of the primary cosmic rays. The actual 


existence of the fine structure has been established 
by experiments of Ribner® and Cooper® in 
Missouri, and of Schremp and Baijios’ in Mexico 
City. 

In Missouri-the fine structure pattern was ob- 
served to have a circular symmetry about the 
zenith. This was interpreted as indicating the 
existence of lines or bands in the primary energy 


* From a dissertation submitted in partial fulfillment of 
the requirements for the Ph.D. degree at Washington 
University. 

t Now at the Radiation Laboratory, Massachusetts 
Institute of Technology, Cambridge, Massachusetts. 

1M. S. Vallarta, rank. Inst. 227, 1 (1939) 
E. J. Schremp, P ys. Rev. 54, 157 (1938); R . Albagli 
Hutner, Ph "have. 55, 15 and 614 (1939). 

*E. J. Schremp, Phys. Rev. 53, 915A (1938). 

. Schremp, Phys. Rev. 54, 157 (1939). 
‘E. . Schremp and H. S. Ribner, Rev. Mod. Phys. 11, 


149 (1 39). 
5H. S. Ribner, Phys. Rev. 56, 1069 (1939). 


: y M. Cooper Rev. 58, 288 (1940). 
614 1941 


‘spectra, for particles too energetic to be ob- 


servably influenced by the earth’s magnetic field.® 
At Mexico City one circular ring of prominences 
was observed, corroborating the results obtained 
in Missouri. Other prominences in the Mexico 
City pattern failed to exhibit this circular sym- 
metry, but were nearly mirror images of one 
another in the N-S plane. The latter provided 
evidence that the primary radiation contains posi- 
tive and negative particles with the same e/m.* 

Although these experiments dealt with the 
static characteristics of the fine structure pattern, 
the existence of time variations in the pattern 
was suspected.'® In spite of the realization that 
new information might be revealed by the study 
of these time variations," the matter could not be 
pursued at that time because of limitations in 
apparatus. The first attempt to explore such 
cosmic-ray fine structure time variations, with 
new apparatus, is discussed in this paper. 


II. APPARATUS 


A directional detecting system with a high 
counting rate combined with high angular resolu- 
tion is needed. The cosmic-ray telescope used in 


- this experiment was designed to measure radia- 


§ See reference 4, p. 150. 
* Cf. reference 7. 

10 Cf. references 4-6. 

1 Cf. references 4-6. 
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tion from eighteen directions simultaneously, 
using fairly large Geiger counters and large 
spacings between counters. The cross section for 
each of the directional systems is about 2.5 
square feet. The dimensions are so adjusted that 
the geometrical angular aperture for any direc- 


Fic. 1. Multidirectional cosmic-ray tel . At the 
time this picture was taken only eighteen of the seventy- 
two Geiger counters were in position on the drum. 


tion is 7° 40’ in the plane of the drum and 18° 18’ 
in the plane of the counter anode wires. 

The use of this equipment, which is described 
in detail below, has indicated how certain im- 
provements in apparatus would enable one to 
secure more complete information from this type 
of experiment. A summary of these suggested 
improvements in apparatus is included in 
Appendix III. 


1. The Geiger-Miiller Counters 


The type of counter used. is essentially that 
recommended by Trost,” consisting of a brass 
cathode, a tungsten anode, and a filling of argon 
plus organic vapor. A departure from conven- 
tional counter construction lies in the elimination 
of the usual glass envelope. The filling gas in- 
cludes a mixture of alcohol and petroleum ether 
which minimizes temperature effects in the 
counting rate. 


® Adolf Trost, Zeits. f. Physik 105, 399 (1937). 
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Further details of the Geiger counters are given 
in Appendix I. 


2. Mechanical Arrangement of Telescope 


The multidirectional cosmic-ray telescope is 
shown in Fig. 1. The Geiger-Miiller counters are 
mounted on the outer periphery of a drum six 
feet in diameter. One Neher-Pickering quenching 
circuit is provided for each set of counters; these 
circuits are mounted on a side panel which circles 
the drum on one side just below the counters. 
Shielded cables, encircling the drum in a layer 
beneath the counters, distribute pulses from the 
quenching circuits to the mixing circuits, which 
are below the cables. The high voltage power 
supplies are in a cabinet built into the base; the 
eighteen message registers which record the im- 
pulses for the eighteen directions are also in this 
cabinet. Slip rings and brushes link the circuits 
on the drum and in the cabinet. In the cabinet all 
the meters are on one panel; a movie camera is 
set up to photograph this panel automatically. 

Changing the azimuth in which the plane of 
the drum lies can only be accomplished by 
moving the whole structure of the telescope. The 
zenith angle, however, can be changed auto- 
matically. The drive for this is provided by a 
$-hp, 1725-r.p.m. G.E. motor, with a 900:1 
Boston “Reductor” and a 60:1 worm and gear. 
The automatic control circuit, for the drum and 
the camera, is in the cabinet with the power 
supplies. 


3. Electric Circuits 
a. Automatic Control Circuit — 


This circuit carries out these functions: It 
automatically starts the motor at periodic time 
intervals of twenty minutes, thirty minutes or 
one hour, and automatically stops the motor 
when the drum has rotated through a prede- 
termined angular interval of one, three, or six 
degrees; it operates a message register which 
indicates the new position of the drum after the 
rotation, and automatically operates a camera to 
photograph the readings of the recording instru- 
ments at the end of the rotation of the drum. 


Provision is made for leaving the drum stationary - 


and automatically photographing the instrument 
panel at intervals of twenty minutes, thirty 
minutes, or one hour. 
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we 


+ 


3... TO FIL. OF 
POWER SUPPLIES 


switcH) 
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CLOCK MOTOR 
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Fic. 2. Control circuit providing for automatic rotation of drum and photographing of meters in the multidirectional 
cosmic-ray telescope. 


The essential features of this circuit, and the 
design of the special switch which makes this 
operation possible, have been described else- 
where ;# a simplified circuit is given in this refer- 
ence. Additional explanatory notes for the com- 
plete control circuit, which appears in Fig. 2, are 
given ini Appendix IV. 


b. Neher-Pickering Quenching Circuit 


This is a conventional circuit; its operation is 
described elsewhere.'* The constants are given in 
Fig. 3. 


c. Mixing and Recording Circuit 


This circuit, shown in Fig. 4, is a combination 
of the Rossi coincidence circuit'* and the Johnson 


+1725 


OHM 


NEGATIVE PULSE 
fom R,*20000 OHMS 


Fic. 3. Neher-Pickering quenching circuit. 


4M. L. Yeater, Rev. Sci. Inst. 14, 146 (1943). 


4 A. L. Hughes, Am. Phys. T. 7, 271 (1939). 
18 B. Rossi, Nature 125, 636 (1930). 


multivibrator circuit.’* Its function is to accept 
pulses from five co-directional trains of counters 
(i.e., from the quenching circuits associated with 
these trains) and cause a message register to 
record a count whenever any one of these trains 
experiences a coincidence of cosmic-ray impulses 
in its members. The operation, briefly, is as 
follows: The 6C6 mixer tubes are normally con- 
ducting; whenever two negative pulses from 
counters on opposite sides of the drum come to 
the grids of any pair of the 6C6’s, say A and B, 
the voltage across R,’ rises sharply and fires neon 
lamp Ne’. This causes the drop across R, to 


Fic. 4. Mixing and recording circuit for one direction; 
Ri=R,' =50,000 ohms, megohm, R;=0.1 megohm; 
Ne, }-watt neon lamps. 


1¢T, H. Johnson, Rev. Sci. Inst. 9, 221 (1938). 
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change so that the 89; a power tube, conducts. 
The duration of the conducting state in the 89 is 
increased by the feedback through C so that the 
message register coil MR in the plate circuit of 
the 89 can operate a mechanical numbering 
device. 

Eighteen such circuits are used, one for each 
direction for which cosmic-ray intensities are 
measured. These are mounted on separate panels 
which encircle the drum; they are visible in 
Fig. 1, and have been identified in the description 


already given for the mechanical arrangement of 


the cosmic-ray telescope. The pulses received by 
these circuits from the quenching circuits must be 
sorted by a system of feeder cables, so that each 
mixer tube receives pulses only from the proper 
quenching circuit. This cable system is described 
in Appendix II. 


d. Power Supplies 


The regulated power supply for the plates of 
the 6C6 tubes in the coincidence circuits follows 
the circuit of Fig. 5. The screens of these tubes are 
provided with a similar supply which is also used 
for the negative voltage for the bias of the 89 
power tube. The third supply, which operates the 
message registers through the power tube in the 
recording circuit, also uses the same circuit; two 
2A3’s are used in parallel, however, and heavier 
transformer and choke are used, because of the 
larger current. 

The circuit shown in Fig. 6 provides the high 
voltage for the Geiger-Miiller counters. In the 
type of experiment here involved, where the 
apparatus must run continuously for weeks or 
months, the first requirement for the high voltage 
supply is that it exhibit no drift in output voltage. 
It was found desirable to take these precautions: 

(1) In order to avoid possible changes in the 
grid voltage for the 6J7, in the regulator circuit, 
the usual potentiometer was not used lest a 
defective contact arise at a, Fig. 6, where a 
changing potential is normally the means of con- 
trolling the output voltage. 

(2) Batteries were used in the cathode circuit 
of the 6J7, in the regulator; neon lamps were first 
tried but were not reliable. 

(3) The vacuum tubes are all mounted in the 
upright position. The voltage becomes noticeably 
erratic if this is not done. 


When these precautions were taken the drift 
amounted to 1.5 percent at most; some of this 
may have been meter drift. No ripple greater 
than one volt could be detected in the output by 
the oscillograph. 

As an additonal precaution, the input to all 
power supplies came from heavy duty Sola con- 
stant voltage transformers. 


4. Light and Temperature Sensitivity 


A serious difficulty arose from the effect of 
light on the sensitivity of the system. A periodic 
daily change in counting rate was observed. This 
effect was traced to the neon lamps in the 
multivibrator recording circuits. Coating them 
with Glyptal lowered the counting rate of the 
system by a factor of ten or more, but without 
changing the ratio of vertical to horizontal count- 
ing rate for one train (one directional pair of 
counting units). The recording circuits were very 


CH 
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Fic. 5. Typical circuit for power supplies used for mixing 
and recording circuit. 
R,=10* ohms, 2 watts. 
R:=25,000 ohms, —2 watts. 
R;=15,000 ohms, wire wound potentiometer. 
R,= 5000 ohms, 1 watt. 
R=50,000 ohms, 1 watt. 
Rs= 250,000 ohms, 1 watt. 
C, =0.002 f, tubular. 
C:=0.01 f, tubular. 
C;=C, 16 f, electrolytic. 
Ch=Thordarson choke, T-74C29. 
T = Thordarson, T-13R08. 
(For power tube plate power supply: Ch-T-74C29, 
T-13R09). 


inefficient with neon lamps in darkness (normally 
they always received a little light at night from 
the light used in conjunction with the camera). 

This difficulty was met by the direct method of 
putting extra lamps in the tent which housed 
the telescope, to provide continuous diffused 
illumination. 

Ambient temperature changes for the cosmic- 
ray telescope were pronounced, as it was housed 
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in a tent on the open roof of the Physics Building 
at Washington University. Temperature readings 
were recorded several times each day. No corre- 
lation between temperature and counting rate 
was discovered for temperatures above 20°F. Any 
temperature effects should affect all directions 
uniformly, and corrections were applied to the 
data to take account of possible aberrations of 
this type. 
Ill. PROCEDURE 
The system was operated continuously, except 
when the temperature was below 20°F. Every 
two to eight hours the drum was rotated by hand 
through two or two-and-one-half degrees; rota- 
_ tions were continuous in the same direction. The 
meter panel was photographed automatically 
every half-hour. Since the directional counting 
arrays are separated in angle by ten degrees, 
_either five two-degree rotations or four two-and- 
one-half degree rotations provided data covering 
the complete 180-degree range in zenith angle. 
This complete coverage was accomplished in 
periods of from twelve to twenty-eight hours. 
The advantages of this manual rotation of the 
telescope drum over automatic rotation are 
several, and a summary of these will serve to 
bring out more information concerning the experi- 
mental method used in this work: (1) The func- 
tioning of the telescope can be checked at the 
time of rotation of the drum, which is when 
unfortunate minor accidents are most likely to 
occur. (2) The routine calculations, exclusive of 
finding probable errors, are reduced by lengthen- 
ing the time intervals during which the drum 
remains in one position ; with automatic rotation 
this would have been one hour at most. (3) The 
total time interval for each zenith direction is 
broken up into subintervals by the automatic 
photographs, so that one can check the probable 
error of any point. (4) The number of counts per 


YEATER 


Fic. 6. High voltage power supply with regulator. 


point on a complete intensity curve is statistically 
more significant than would be the case if the 
drum were at one position for only one hour or 
less. All of these advantages except (1) could be 
achieved with completely automatic operation of 
this apparatus if the control circuit were modified 
to permit the use of other angular intervals and 
other time intervals. 


Whether operated manually or automatically, — 


the present type of apparatus allows the angular 
interval between successive positions of the drum 
to be so chosen that each directional array is 
made to explore one. ‘‘section’’’ of the zenith 
angle range. These sections can be joined to form 
one continuous curve, as described below. 

Readings of temperature and atmospheric 
pressure were recorded several times daily while 
cosmie-ray observations were carried out. 

The data were taken from the photographic 
record as soon as possible and later reduced as 
described below. 

A few improvements in procedure which this 
work showed to be desirable are suggested in 
Appendix III. 


IV. TREATMENT OF EXPERIMENTAL DATA 
1. Corrections 


It was considered advisable to make certain 
corrections. Since the western end of each section 
in the east-west curves, for example, begins at 
the same time, a changing counting rate caused 
by spurious effects common to the whole instru- 
ment (such as voltage changes or ambient 
temperature changes) could lead to spurious 
anomalies common to all sections. A good meas- 


-ure of the drift in counting rate is afforded by the 


total counting rate for all eighteen trains. In the 


17 The word “section” will henceforth, in this paper, be 
used in this sense, denoting a small portion of 
angle range. 
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absence of spurious effects this total counting 
_ rate should be nearly the same for every position 
traversed by the drum within a ten-degree sec- 
tion. The method of making the corrections, 
then, involves multiplying each recorded number 
of counts by such a factor that the total corrected 
number of counts for each position of the drum is 
the same and equal to the average total value 
calculated from the observed totals. The multi- 
plying factors are just the ratios of the various 
observed total counting rates, one for each posi- 
tion of the drum, to their average value for the 
given ten degree section. The same factor applies 
to every one of the eighteen intensities recorded 
for each position of the drum. 

These results of this method of applying cor- 
rections are to be noted: (1) Any common 
bending of the ten-degree sections, which might 
lead to spurious anomalies, will be straightened 
in this process. (2) Any directional effects of 
atmospheric or magnetic disturbances will show 
up only in the directions actually involved. For 
the anomalous change in one direction will be a 
small fraction of the total intensity in that direc- 
tion, which in turn will certainly not be a large 
fraction of the total intensity for eighteen direc- 
tions. (3) These corrections do not eliminate 
errors in relative intensity between any two 
sections, caused by different counting rates for 
different directional arrays. Although the sections 
can be adjusted to form one continuous curve, as 
described below, errors can accumulate in this 
‘ curve in the process of bringing the separate ten- 
degree sections into coincidence at their terminal 
points. For a statistical error, say, in one terminal 
point may be introduced proportionally into 
every point of the adjacent section. This means 
that different sections, either before or after being 
merged into one continuous curve, can be com- 
pared with respect to shape but not with respect 
to intensity; and in particular, that east-west 
asymmetries cannot be determined from the 
present data. 


2. Adjustment of Sections to Form 
Complete Curves 


The system was rotated through four 2.5° 
intervals or five 2° intervals to secure data for one 
curve. Since the individual arrays are separated 


in angle by ten degrees, sections thus found have 
common end points. The adjustment of the 
sections to form a continuous curve involves 
multiplying each point of a section by that factor 
which makes one of its terminal points coincide 
with the terminal point of one adjacent section. 
The other adjacent section is then similarly 
matched to this section. 


3. Probable Errors 


Since automatic recording of meter readings 
occurred every half-hour, each time interval for 
one angle is effectively broken into subintervals. 
From these data one could compute the probable 
error of every point by usual methods.'® Since an 
examination of the data indicated no excessive 
fluctuations, the labor of such calculations has 
been omitted. An estimate of the probable error 
of any point is given by m', where n is the number 
of counts. The statistical treatment described 
below, moreover, relieves the probable errors of 
final responsibility for the validity of the results. 


4. Statistical Analysis 


In order to find time variations in the direc- 
tional intensity pattern, the daily curves are 
compared with the time-average of all such curves 
for the same azimuth. Pearson’s criterion for 
goodness of fit'® is used. Each ten-degree section 
of an individual curve is compared with the 
corresponding section of the time-average curve 
after the latter has been fitted as follows: Let me 
be the number of counts per point, at the zenith 
angle 0, of the time-average curve. Let rp be the 
number of counts per point, at angle 6, of the 
individual curve. Then 


re=ne 


(0;=0°, 10°, 20°---, or 2°, 12°, 22°---, etc.) 


gives the points for the time average section 
fitted to the individual section. The deviation 5» 


18 Cf. Gunnar Dah , Statistical Methods for Medical 
and Biological Students ( Allen and Unwin, London, 
1940), Chap. VIII; T. C. Fry, Probability and Its Engineer- 
rr cGraw-Hill Book Company, Inc., New York, 

19 Discussed by Fry, reference 18, p. 285. 
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for any point of the individual curve is then 
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Fic. 7. Time average cosmic-ray directional intensi 
I(@) as a function of zenith angle @ for the E-W azimu 
(upper curve) and N-S azimuth (lower curve). For the 
E-W curve each vertical scale division represents 2000 
counts; for the N-S curve each division represents 1000 
counts. The maximum intensity for each curve is shown. 
The time interval covered by the E-W curve extends from 
9:10 p.m., March 25, 1943 to 8:10 p.m., April 12, 1943; the 
time interval for the N-S curve is from 3:35 p.m., April 15, 
1943 to 1:06 p.m., April 20, 1943. The average probable 
error is about 0.4 percent of the maximum intensity for the 
E-W curve and 0.7 percent of the maximum intensity for 


the N-S curve. 


one enters Pearson’s table for goodness of fit 
criteria.” The probability number P given by the 

2° Karl Pearson, Tables for Statisticians and Biometricians 
(Cambridge University Press, London, 1930), third edi- 
tion, Part I, p. 26. 
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table is the probability that another run of the 
experiment would by chance yield a curve 
deviating from the standard of comparison as 
much or more than the curve at hand. In the 
present data each ten-degree section of a curve 
contains either five or six points; one therefore 
enters Pearson’s table with n’° equal to four or 
five, respectively, since a relationship between 
the points is implied by the method of fitting the 
time-average curve to the individual curve. 


= 
60 70 60 50 40 30 20 10 O 10 20 30 40 50 60 70 80 


Fic. 8. Individual short-time directional intensity curves 
for the north-south azimuth. Each vertical scale unit 
represents 250 counts. The maximum intensity is shown for 
each curve. The average probable error is about 2 percent 
of the maximum intensity. The time intervals are as follows 
(all dates for 1943): I. 4-15, 3:35 p.m. to 4-16, 6:06 a.m. 
II, 4-16, 9:06 a.m. to 4-16, 7:06 p.m. III. 4-16, 11:36 P.M. 
to 4-17, 3:06 p.m. IV. 4-17, 3:36 P.M. to 4-18, 1:36 A.M. 
V. 4-18, 3:36 a.m. to 4-18, 7:36 p.m. VI. 4-18, 8:36 P.M. to 
4-19, 8:36 a.m. VII. 4-19, 4:36 P.M. to 4:20, 2:36 a.m. 
VIII. 4-20, 3:06 a.m. to 4-20, 1:06 P.m. 


V. RESULTS AND DISCUSSION 


1. The Time-Average Directional 
Intensity Curves 
For the east-west and north-south azimuths 
the cosmic-ray intensity J(@) is plotted as a 
function of zenith angle @ in Figs. 7-9. The two 
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curves of Fig. 7 are summation, or time-average, 
curves representing all of the data recorded for 
the E-W and N-S azimuths, respectively; they 
are obtained by the superposition of the short- 
time curves of Figs. 8 and 9, which will be dis- 
cussed in succeeding sections. 

These time-average curves are similar in scope 
to the measurements made previously by Ribner™ 


COSMIC-RAY INTENSITY 


Fics. 9A-C, Individual directional intensity curves for 
the east-west azimuth. Each vertical scale unit represents 
500 counts for curves I-IV and XIII-XVI, and 250 counts 
for curves V-XII. The maximum intensity for each curve is 
shown. The average probable error is about 1 nt of the 
maximum intensity for curves I-IV and XIII-XVI, and 
about 2 t of the maximum intensity for curves 
V-XII. The time intervals are as follows van dates for 
1943): I. 3-25, 9:10 p.m. to 3-27, 1:10 a.m. II. 3-27, 2:40 
A.M. to 3-28, 6:40 a.m. III. 3-29, 4:40 p.m. to 3-30, 7:10 
p.M. IV. 3-30, 7:10 p.m. to 3-31, 9:10 p.m. Each of curves 
V, VII, [X, and XI covers 12 hours starting at 7 P.M. on 
3-31, 4-1, 4-2, 4-3, respectively. Curves VI, VIII, X, and 
XII are for 12-hour periods starting at 7 A.M. on 4-1, 4-2, 
4-3, 44, respectively. XIII. 4-7, 3:10 a.m. to 4-8, 8:10 
AM. XIV. 4-8, 11:40 a.m. to 4-10, 12:10 a.m. XV. 4-10, 
Heo A.M. to 4-11, 12:10 p.m. XVI. 4-11, 5:10 p.m. to 4-12, 

P.M. 


and Cooper” in Missouri, although more cosmic- 
ray counts are recorded here than in the earlier 
work. This permits a comparison of the present 
results with those previously reported. 

In the present time-average curves prominences 
can be recognized near E 5°, E 10°, E 20°, W 5°, 
W 10°, W 20°, S 5°, S 10°, S 20°, and N 5°. At 
N 20°, where others have reported an anomaly, a 
prominence is just perceptible. Prominences at 
40° in the north, south, east, and west are not 
strong but are sufficiently evident to bear out the 
earlier reports that anomalies do exist at these 
angles. In addition to the 40° anomaly, the data 
of Ribner and Cooper reveal prominences near 7° 
and 20° in the E and W azimuths (Ribner and 
Cooper) and in the N, S, NE, NW, SE, and SW 
azimuths (Cooper). The prominences in the 


*1 See reference 5. 
% See reference 6. 


the 
irve —3700 _ 
1 as | 116) 
the _3500__ xu 
irve 1 
fore 
r or 
Il 
reen 
the 
Vv 
w 
680 70 60 SO 40 30 20 10 © 10 20 30 40 50 60 70 8° 
Fic. 9A, 
Fic. 9C. 
1500 
Vv 
vi 
1300 _ 
1750 
1x 
j _!550__ 
BO 
x! 
irves 
n for 
cent 
lows / 
A.M. 
P.M. 1 
A.M. 
{. to 
AM. 
iths 
sa 


82 MAX L. YEATER 


present curves near 5° and 10° are probably the 
same as appeared to be one near 7° from the more 
widely spaced points of the earlier work. The 
prominences at 20° in the N, S, E, and W appear 
to be the same for both the present and past data. 

The present time-average curves, therefore, 
tend to confirm rather completely the existence 
and circular symmetry of the anomalies previ- 
ously reported; they support the interpretation 
suggested at that time for this circular symmetry, 
namely, that lines or bands exist in the primary 
energy spectrum.” 


2. The Daily Curves and Their Time Variations 


Each of the twenty-four curves of Figs. 8 and 9 
spans a period of from twelve to twenty-eight 
hours, and covers the full 180° of one azimuth at 
2° or 2.5° intervals. An attempt was made to 
determine the existence and character of time 
variations in the directional cosmic-ray intensity 
merely from an inspection of these daily curves. 
This rather obvious step was found inadequate to 
provide any consistent interpretation of the data; 
selected combinations of these curves were 
similarly uninformative. 

A better approach to the study of these curves 
is a section by section comparison of corre- 
sponding ten-degree angular intervals of the 
different curves ;* this is best accomplished by a 
statistical analysis. The method adopted has 

3 Cf. references 4, 6, and 7. 
*% For this purpose the original corrected data for the 


individual ten-degree sections are used, rather than the 
adjusted. values which form the continuous curves of 


Figs. 8 and 9 


already been described in Section IV, 4, and the 
results of this analysis will now be discussed. 
For the east-west azimuth, Table I gives the 
probability values P describing the goodness of 
fit of each ten-degree section of the individual 


. directional intensity curves with the time-average 


curve; Table II lists the corresponding values for 
the north-south azimuth. There are many entries 
for P<0.1 and, in the E-W azimuth, even for 
P<0.01. In determining whether these low values 
actually indicate significant time variations, three 
interrelated questions are to be dealt with: (1) A 
few low P values are to be expected just by 
chance; is the number of these values here so 
high as to imply the existence of real time 
variations? (2) What limiting value of P shall be 
considered useful for distinguishing abnormal 
sections of the individual curves, and hence for 
indicating the angular distribution of any moving 
prominences? (3) Does the angular distribution 
of the low P values indicate the absence of 
random instrumental errors? 

The curves of Figs. 10 to 16 are helpful for a 
consideration of these questions, which will be 
taken up in order. 

(1) In each case of Figs. 10 to 13, for the E-W 
azimuth, the horizontal ‘chance line’ is ap- 
preciably below the actual distribution. This line 
becomes relatively lower, furthermore, as the 
limiting values of P are decreased from 0.2 to 
0.01; the average ratios of chance expectation to 
actual number of sections are 0.6 for P<0.2, 0.4 
for P<0.1, 0.3 for P<0.05 and 0.17 for P<0.01. 
The existence of significant time variations in the 


TaBLE I. Probability values P describing the goodness of fit of the individual east-west curves with their time average.* 


Ill IV v VI 


XII XIII XIV XV 


x XI XVI 


geek 


0.50 0.70 
16 13 


* The missing entries for this table represent data thrown out because of instrumental errors. 
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NUMBER OF SECTIONS 


85 75 6S 55 45 35 25 15 5 S '5 25 35 45 55 6S 75 85 


w ZENITH ANGLE , DEGREES e 
Fic. 10. Number of 10° sections of individual E-W directional intensity curves for which P=0.2, 
as a function of zenith angle. The broken horizontal ‘“‘chance line” shows the most probable number 
expected by chance for any 10° interval. 


- 


NUMBER OF SECTIONS 


3218S 25S 3 3S GS 7S SS 
w ZENITH ANGLE, DEGREES 
Fic. 11. Number of 10° sections of individual E-W directional intensity curves for which P=0.1, 


as a function of zenith angle. The broken line shows the most probable number expected by chance 
for.any 10° interval. 


east-west directional intensity pattern therefore instrumental errors can be shown to be absent. 
seems to be indicated, provided only that random The results for the north-south azimuth are 
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NUMBER OF SECTIONS 


85 75 6S 55 45 35 


235 15 3 5 


iS 25 35 45 55 6 75 85 


Ww ZENITH ANGLE , DEGREES Ee 


Fic. 12. Number of 10° sections of individual E-W directional intensity curves for which P=0.05, 
as a function of zenith angle. The broken line shows the most probable number expected by chance 


for any 10° interval. 
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Fic. 13. Number of 10° sections of individual E—-W directional intensity curves for which P=0.01, 
as a function of zenith angle. The broken line shows the most probable number expected by chance 


for any 10° interval. 


different. For P<0.2 and P<0.1 the number of 
sections is in no case appreciably greater than 
could be expected by chance (cf. Figs. 15 and 16). 
Only two sections exhibit a value less than 0.05, 
although seven would be the normal chance 
expectation; none has P<0.01. This situation 
will be discussed further after questions (2) and 
(3) have been taken up. For the moment we 
merely point out this difference between the data 


for the N-S and E-W azimuths, a difference 


which implies that time variations of a type which 
can be revealed by Pearson’s test exist in the 
E-W azimuth but not in the N-S azimuth. 

(2) It is reasonable to use the plot of P<0.05 


as a guide in locating the angular distribution of 
time-varying sections of the pattern. This choice 
can be justified to some extent by the ratio of the 
number of sections expected by chance to have 
P<0.05, to the number actually observed, for the 
E-W azimuth (cf. Fig. 12). The average value of 
this ratio is 0.3 and even lower values exist for 
several angles, the minimum being 0.11. The 
corresponding ratio (average value equals 0.17) 
for P<0.01 of course does not preclude 0.01 as 
an upper limit of P, but the observed number of 
sections for this case is too small to be useful. 
(3) The angular distribution for P<0.05 in 
the E-W azimuth suggests that instrumental 
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errors are relatively insignificant. For low P 
numbers caused by instrumental aberrations may 
be expected to show a rather uniform angular 
distribution. This is not the case in the present 
data, as shown by Fig. 12. Of particular im- 
portance are the 80°-90° sections, where the 
“accidental’’ counting rate of the telescope should 
contribute the greatest relative intensity. Of 
thirty-two sections for this range, only one ex- 
hibits a P<0.05; none a P<0.01. That instru- 
mental errors are unimportant is further indicated 
by the absence in the north-south data of many 
low P numbers. 

We return now to the discussion of the differ- 
ences between the N-S and E-W observations. 
It is important to know if this difference is a real 
cosmic-ray effect or if it can be accounted for by 
limitations in the experimental or analytical 
procedure.*® These possibilities must be con- 


% A re-examination of the data by another statistical 
method is desirable, but up to the date of this writing there 
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sidered : (1) That the total time interval covered 
by the N-S data did not include magnetic or 
barometric changes like those which occurred 
while the E—W observations were made. (2) That, 
in view of the smaller number of curves for the 
N-S azimuth, there should be effectively a 
fractional increase in the number of constraints*® 
pertaining to these data in the application of 
Pearson’s test; for an implicit relationship be- 
tween the reference (normalized time-average) 
curve and any individual short-time curve exists, 
since the reference curve is an average of all of the 
individual curves. (3) That a real difference in 


has been no time for this. In a forthcoming paper by 
Schremp and Yeater, however, reference wil be made to 
the results of this second analysis, which will have been 
carried out by then. 

26 The number of constraints for a section in either 
azimuth has been taken as one; this constraint enters be- 
cause the total number of counts in the reference section is 
made equal to the total number in the section being tested 
for a time variation. ‘“‘Constraint,” in the sense implied 
here, is defined in Fry, reference 18. 
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. Time distribution of the 10° sections, of individual E-W curves, for which P=0.01. The 
time of day indicated for each date is noon. 
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Fic. 15. Number of 10° sections of individual N-S directional intensity curves for which P=0.2, 
as a function of zenith angle. The broken horizontal “chance line” shows the most probablejnumber 


expected by chance for any 10° interval 
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the nature of the time variations in the cosmic- 
ray directional intensity distribution exists be- 
tween the N-S and E-W azimuths. All of these 
possibilities will now be examined. 

(1) The complete time interval for the eight 
N-S curves includes five days as against eighteen 
days, March 25 to April 12, 1943, for the sixteen 
E-W curves. During these five days, from April 
15 to April 20, 1943, the international magnetic 
character numbers (cf. Fig. 17) show no ex- 
ceptional uniformity; no general difference be- 
tween the magnetic fluctuations associated with 
the N-S and E-W cosmic-ray data is revealed. 
There is the particular difference, however, that 
no magnetic disturbances severe enough to be 
classed as storms occurred between April 15 and 
April 20; while three brief storms between March 
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23 and April 12 were reported by some of the 
stations making magnetic observations.”” As to 
the atmospheric pressure, the barometer readings 
at St. Louis for the period April 15 to April 20 
show a somewhat more stable condition than 
existed at several times during the preceding 
month (cf. Fig. 17). A change of pressure 
amounting to 0.6 cm Hg did occur, however, 
while the N-S cosmic-ray data were taken; the 
changes were not more than 1.0 cm throughout 
most of the time for the E-W data. 

The possibility exists, therefore, that the N-S 
results may be caused by the somewhat more 
quiescent conditions of the magnetic field and/or 
atmospheric pressure ; but the evidence favoring 
this possibility is not entirely convincing. For the 


*7 Terr. Mag. 48, 117, 188 (1943). 
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Fic. 16. Number of 10° sections of individual N-S directional intensity curves for which P=0.1, ~ 
as a function of zenith angle. The broken line shows the most probable number expected by chance 
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' Fic. 17. Variation of atmospheric pressure at St. Louis and of international magr 
number, for the time covered by the present cosmic-ray data. The time of day indicated for each 
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difference between the N-S and E-W cases, with 
regard to magnetic and barometric conditions, is 
not great. 

(2) Even if we take the somewhat extreme 
steps of increasing the number of constraints for 
the N-S case only, and letting this increase equal 
one, the difference between the N-S and E-W 


results is not removed. In this case an entry of 


P=z0.1 in Table II becomes P <0.05 ; Fig. 16 then 
becomes a graph of P<0.05, with the height of 
the ‘chance line’’ reduced from 0.8 to 0.4. For 
this assumed condition the average ratio of 
chance expectation to observed value for the 
number of sections with P< 0.05 is 0.72, and also 
for P<0.01 it is 0.72, with no entries for P<0.01; 
these are to be compared to values of 0.3 and 
0.17, respectively, for the corresponding ratios in 
the E-W case, where in addition P values lower 
than 0.0001 are found (Table 1). 

From this it appears that the difference in the 
present results for the N-S and E-W azimuths 
_ has its origin in some physical circumstance. 

(3) This possibility, that the character of the 
time variations in the cosmic-ray directional in- 
tensity patterns is a function of azimuth, is not 
ruled out by the preceding considerations (1) and 
(2). That time variations in the N-S pattern may 
be absent is consistent, furthermore, with certain 
interpretations of the E-W data, as discussed in 
the next section. This third possibility, in the 
light of all the evidence, appears by far the 
strongest. 


3. Implications of the Time Variations 


One of the purposes of this experiment was to 
determine whether rapid time variations exist for 
the anomalies observed in the static cosmic-ray 
directional intensity pattern. The result of the 
experiment in this respect is that no significant 
variations in these anomalies are revealed by the 
present analysis. On the other hand, this analysis 
does indicate time variations of an independent 
character. 

The immediate evidence for the independence 
of the static and dynamic anomalies lies in the 
fact that the circular symmetry exhibited by the 
static anomalies is not exhibited by the time 
variations in the pattern. The fundamental 
reason for this independence, however, would 
appear to be that different origins can be ascribed 


to the two effects. The circularly symmetrical 
anomalies in the time-average curves have their 
origin in atmospheric effects, as previously dis- 
cussed by others.** The probable origin of the 
time variations will be discussed in the following 


paragraphs. 
In the E-W azimuth a dynamic character in 


TaBLE II. Probability values P describing the go 
eir time 


of fit of the individual north-south curves with 
average. 
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the cosmic-ray directional intensity pattern is 
indicated by Fig. 12 to exist near W 65°, and 
E 75°, and near W 15° and E 25°. As we shall see, 
there is reason to believe that these peaks actu- 
ally represent symmetrical dynamic regions. For 
if the variations in the pattern which produce the 
peaks in Fig. 12 really occur at 20° and 70°, re- 
spectively, then the observed asymmetry can 
easily result from a slight inaccuracy in the de- 
termination of the vertical direction for the 
cosmic ray telescope. An error of even a degree or 
two in this determination®® can cause a time 
variation at 20° to appear at 25° on one side and 
15° on the other in Fig. 12, and similarly for 
variations at 70°; for a point at E 25°, for ex- 
ample, means only that the variation in the © 
pattern lies between 20° and 30°. Such an 
explaining away of the small observed asymmetry 
of course requires that this asymmetry be in the 


8 Cf. references 4, 6, and 7. 
** The geomagnetic vertical should be used when mag- 
e geographic vertical at St. Louis; the geographic 
vertical, however, is used for the present data. If the 
pega vertical is used, the s in Fig. 11 tend to 
me symmetrical. The experimental alignment of the 
— . with the vertical direction was accurate to 
a 


he 
to 
igs 
20 
an } 
ng | 
ire 
er, 
he 
ut 
Curve I u oom Vv VI Vl 

Angular 
re section | 
| S 90-80 
= } 
ng 60-50 
he 
30-20 
20-10 | 
10-0 

N 0-10 

10-20 
20-30 
30-40 

40-50 

50-60 

60-70 
70-80 
80-90 | 
| 
| 


88 


same direction for both the 20° and 70° cases; 
this requirement is fulfilled. 

Another feature to be noted in Fig. 12 is that’ 
the higher peaks occur at the larger angles. The 
greater tendency to change thus exhibited at the 
large zenith angle of 70° is an indication that 
magnetic effects, rather than atmospheric pres- 
sure changes, are responsible for these time 
variations. This follows because variations caused 
by atmospheric pressure changes should be more 
pronounced at smaller angles.*® Variations caused 
by magnetic changes, as brought out pet, ae 
not restricted in this way. 

Further evidence that these cosmic-ray time 
variations are associated with terrestrial magnetic 
changes comes from the time distribution of the 
very low values of the probability index P 
(Fig. 14). It happens that maxima in this curve, 
centered at March 26, April 1, and April 10, are 
observed coincident with or just following the 
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Fic. 18. Simple shadow cones, for positive particles, at 
geomagnetic latitude north 50°. Energies are in Stérmers. 


%° This can be demonstrated as follows: Let he be the 
vertical height of the atmosphere, and h the atmospheric 
thickness looking in the direction @ from the zenith. Then 


to a good approximation 
h=hosec@ or cos0=ho/h. 
Hence 
= — (1/h) cosec — (1/ho) cotan 
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occurrence of magnetic disturbances. (Magnetic 
data are plotted in Fig. 17; the disturbances 
referred to are listed specifically elsewhere.) 
These magnetic changes can affect the direc- 
tional pattern by distorting the allowed cones of 
the primary cosmic rays. In considering the effect 
of changes in these cones we perceive that the 


absence of time variations in the N-S pattern can 


be consistent with the hypothesis that variations 
at other azimuths are caused by magnetic 
fluctuations. 

For directional cosmic-ray changes of magnetic 
origin must be caused by the movement of an 
allowed cone boundary across the direction of 
observation. Now at the magnetic latitude of St. 
Louis, approximately N 48°, the boundary of the 
allowed cone is essentially the common boundary 
of the simple shadow cone and the Stérmer cone. 
The character of this boundary for various 
energies is shown in Fig. 18.** From this we can 
see how certain of these boundaries can produce 
the observed dynamic regions in the E—-W pat- 


tern (at approximately +20° and +70°) without 


intersecting the N-S meridian. At W 20°, for 
example, lies a cone boundary for a positive 
particle with an energy slightly greater than 0.2 
Stérmer. The N-S meridian is in the forbidden 
region for this energy ; hence a movement of this 
cone can produce no changes in the N-S pattern. 
The region near E 20° can be similarly described 
with reference to a negative primary particle. In 
dealing with the region near E 70° we can invoke 
the cone for a positive particle with energy of 
about 0.6 Stérmer. In this case the north-south 
meridian is in the region of full light, and no 
variations in the N-S pattern are to be expected 
from changes in this cone. The corresponding 
cone for a negative particle can account in a 
similar way for the dynamic character of the 
pattern near W 70°. 

In thus relating the observed time variations 
with magnetic fluctuations, we have a reasonable 
explanation of their symmetry in the E-W 
azimuth and their absence in the N-S azimuth. 
More than this, we find in the symmetry of these 
variations a new kind of evidence supporting a 
conclusion previously advanced by Schremp,* to 


3t See reference 27. 
% Cf. E. J. Schremp, toon Rev. 54, 169 (1938). . 


' % Reference 32, Fig. 5 
Cf, reference 
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the effect that positive and negative particles 
_ with the same e/m exist among the primary 
cosmic rays. 
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APPENDIX I 
Geiger Counter Construction 


The cathode is shown in Fig. 19; this is entirely of brass. 
The end fittings marked f are used in the vacuum seal. The 


eum PYREX 
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Fic. 19. Structure of the Geiger counters. A, cathode; 
B, vacuum seal; C, glass lead-ins; D, anode assembly. 
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method of securing this seal is shown in detail in Fig. 19B; 
R is a gum rubber washer which is compressed against the 
inner surface of f and around the glass lead-in g to form the 
seal; the nut N bears on the metal washer W which crushes 
the rubber. The cathode is formed by hard-soldering the 
threaded fitting f to the disk d, then soft-soldering d to the 
light brass cylinder c. The glass lead-in is shown enlarged 
in Fig. 19C. The complete anode assembly is shown in 
Fig. 19D. 

The anode is prepared as follows: A length of 0.010” 
tungsten wire is cut off and a loop formed at each end by 
twisting under a tension of about two kilograms; the wire 
is then rubbed vigorously with fine sand paper. It was not 
found necessary to flash the anodes in vacuum. 

The cathodes are cleaned in an acid solution consisting of 
about three parts HNO; (sp. gr. 1.403), one part H:SO, 
(sp. gr. 1.84) and twenty parts water. After being washed 
several times in tap water and in distilled water the 
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Fic. 20. A, Geiger counter characteristic curve; B, counting 
rate of Geiger counter at different temperatures. 


cathodes are permitted to stand in air a few days to form 
a thin oxide layer on the brass. 

The filling gas mixture consists of 8 mm alcohol, 4 mm 
petroleum ether, and 7.8 cm argon. A typical characteristic 
curve is shown in Fig. 20A. The effect of temperature on the 
counting rate of an early test counter is shown in Fig. 20B. — 
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APPENDIX II 
Pulse Cable Connection 


The manner in which the counter pulses are properly 
distributed to the mixing circuits is indicated in Fig. 21 and 
Fig. 22. Figure 21 is a simplified diagram showing the 
interconnections for just two of the eighteen directions. 
From each quenching unit a shielded cable runs one-fourth 
of the way around the drum on either side, and from these 
cables the pulses are tapped off to the coincidence circuits. 
The mixing circuits for the direction indicated by the 
broken lines are represented by the ten small dots at the 
left, while for the direction defined by the full lines the 
mixing circuits are at the right. Counter units A and A’ 
feed into vacuum tubes a; and a;’; similarly B and B’ feed 
into 5; and by’, etc. In order to interconnect thirty-six 
counter units for eighteen directions this system is simply 
extended; one cable from each counter unit runs halfway 
around the drum and pulses are tapped off at the appro- 
priate mixing-recording panel. The complete wiring dia- 
gram for this is given in Fig. 22. 


APPENDIX III 


Suggested Improvements In Apparatus and Pro- 
cedure for Cosmic-Ray Directional 
Intensity Measurements 


(1) The general circuit arrangement of the multidirec- 
tional telescope is satisfactory, but a recording circuit not 
using neon lamps should be employed. 

(2) The “fast” type of Geiger counter containing organic 
vapor is not good for this work because of limited life. 

(3) The cosmic-ray telescope should be operated in a 
shack in which constant conditions of temperature and 
illumination can be maintained. 


(4) An output pulse length, for the Geiger counter- 
quenching circuit combination, of the order of a microsecond 
should be achieved, in order to reduce the “accidental” 
counting rate. This, it is believed, can be done by adjust- 
ment of constants in the quenching circuit; a resulting ex- 
cessive decrease in pulse size can be compensated by another 
amplifier stage. This seems more practicable than a triple 
coincidence counting system, for the present type of 
apparatus. 

(5) Observations should be made in rapid alternation in 
the N-S and E-W azimuths (or any other pair at right 
angles to each other) so that in effect patterns for both 
azimuths are recorded simultaneously. 

(6) If measurements like those described in this paper 
are made, a second telescope operating at fixed zenith 
angles should be run concurrently, giving direct readings 
of intensity versus time for various angles. If possible, 
ionization chamber measurements of total intensity should 
be carried along at the same time. 

(7) A high angular resolution, of three degrees or less, 
with respect to zenith angle, is desirable for a cosmic-ray 
telescope used in tracing the fine structure; for a more 
precise knowledge of this pattern, which would immeas- 
urably facilitate the interpretation of the data, appears to 
be easily attainable by this means. A second telescope used 
for direct observation of time variations at fixed angles 
would preferably have a resolution in zenith angle equal to 
the angular spacing between the directional trains, in the 
present case ten degrees. 


APPENDIX IV 
_ Notes on Control Circuit 


These notes pertain to Fig. 2. In reference 13 is given a 
description of the manner of starting the motor M (which 
causes the drum to rotate) at optional time intervals of 


Fic. 21. Pulse cable connections for two directions of observation. The dots 
marked with lower case letters represent mixing tubes in the recording circuit. 
Circles marked with capital letters — Geiger counter units. The angular 
separation of these two directions is 1 
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Fic. 22. Complete pulse cable connections for the multidirectional cosmic-ray telescope. 
The Geiger units are numbered / to 36 around the inner circle. The cables are shown as 
concentric semicircles. Every group of dots represents one panel of mixing tubes. Each group 
is assigned the number of one of the Geiger units of the center train for that direction; for 
example, group 1 represents the mixing-recording circuit for the complete directional array 
of counters which includes the following units in double-coincidence trains: 3 and 17, 2 and 


18, 1 and 19, 36 and 20, 35 and 21. 


twenty minutes, thirty minutes, or one hour, and of the way 
in which the motor is stopped after it has experienced a 
predetermined number of revolutions. This number, and 
hence the angle through which the drum is rotated, is 
chosen by selector switch Sw. At the position 2 of Swe, 
relay 2 will operate, and hence the motor will stop, every 
time either pair of contacts a or d of Sws is shorted by the 
rolling contact* which is driven by M through reducing 
gears. This gear ratio is so adjusted that for a six degree 
rotation of the drum the rolling contact passes once over 
every set of contacts a to f, and these sets are spaced sixty 
degrees apart. Hence for the setting 2 of Sw: the motor is 
stopped after every three-degree rotation of the drum. 
Similarly, for setting 1 the motor is stopped whenever any 
set of contacts a to f is shorted, or for every one degree 


* Cf. reference 13. 


rotation of the drum. Setting 3 of Sw. stops the rotation of 
the drum after six degrees. 

When relay 2 and simultaneously relay 4 are activated, 
thereby stopping the motor, the coils MR and L are 
energized. MR is the coil of a message register which counts 
the number of small rotations of the drum and hence indi- 
cates the angular position of the drum at any time during 


‘a run of the experiment. L is a solenoid which trips a 


camera and causes a photographic record of the meter 
panel to be made. ; 

The camera can be tripped automatically by the clock 
switch without any rotation of the drum. This occurs if Sw. 
is at position 2. 

Switch Sw; is a push-button switch by means of which 
the clock switch can be shorted by the operator, thereby 
enabling him to start the motor M or, if Sw, is at 2, to take 
a picture of the meter panel, at any arbitrary time. 
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The x-ray and electron excitation of silver and cadmium has been studied in the region from 
their thresholds to 3.3 Mev. The thick target x-ray excitation curve for silver shows the thresh- 
old to be 1.18+0.03 Mev and higher activation levels are found at 1.59, 1.95, 2.32, 2.76, and 
3.13 Mev, respectively. The half-life of the metastable state of silver is 40.4+0.2 seconds. Both 
the thick target x-ray excitation curve and the electron excitation curve for cadmium indicate 
that the threshold for these processes is 1.25+0.03 Mev while other activation levels are found 
at 1.68, 2.08, and 2.56 Mev. The energy of the metastable level of cadmium, as determined 
from data on the absorption in aluminum of the conversion electrons emitted in the process, has 
been found to be 195 kev and the half-life of the metastable level is 48.7+0.3 minutes. 


INTRODUCTION 


UCLEAR isomerism was first observed 
when it was found that three radioactive 
periods could be produced by the irradiation of 
bromine with slow neutrons. Since there are only 
two stable isotopes of bromine, namely, 3sBr7® 
and 3sBr*, it was necessary to assign two of these 
periods to a single isotope. 

It was pointed out by von Weizsacker® that 
such a condition may be explained by the as- 
sumption that there is a large difference of spin 
between the two states and hence the state of 
higher energy may exhibit a relatively long life- 
time against y-emission. In 1938, Goldhaber, 
Hill, and Szilard* found that the stable isotope 
In"™* possessed such a metastable level which 
could be excited by fast neutrons and which 
decayed to the ground level with a period of 4.1 
hours. Since this time, eleven stable nuclei 
possessing metastable states have been reported.‘ 

The presence of a low lying metastable state in 
a stable nucleus affords a very accurate method 
for determining a series of the energy levels of 
that nucleus. Inasmuch as the transition from 


the metastable to the ground level is forbidden, . 


because of the large difference in angular mo- 
mentum, the direct transition from the ground 
state to the metastable state is just as strongly 
forbidden. Therefore, to reach the metastable 


1C. H. Johnson and F. T. Hamblin, Nature 138, 504 


(1936); A. H. Snell, Phys. Rev. 52, 1007 (1937). 

2C. F. von Weizsacker, Naturwiss. 24, 813 (1936). 

3M. Goldhaber, R. D. Hill, and L. Szilard, Phys. Rev. 
55, 47 (1939). 

‘For a recent reference list see: G. T. Seaborg, Rev. 
Mod. Phys. 16, 1 (1944) 
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state, the nucleus must first be raised to a higher 
activated level, the spin of which is intermediate 
between that of the ground and metastable 
states. The nucleus can then leave this activated 
level either by direct transition to the ground 
state or it may fall to the metastable level with 
a further change of spin, /,,—/,, and the emission 
of energy, hva—hvm. After the nucleus is in the 
metastable state it can only fall to the ground 
state with the emission of energy hv,» and a large 
change of angular momentum, /, 

Nuclei, in general, can be raised to the activated 
levels by inelastic collision with neutrons,’ pro- 
tons,® alpha-particles,® electrons,’? and by the 
absorption of x-rays. The latter two processes 
afford more complete and accurate methods for 
studying the excitation and for the determining 
of the positions of the higher activation levels. 

If x-rays are produced by a beam of mono- 


chromatic electrons striking a thin target, the . 


intensity of a given isochromat, in the continuous 
x-ray spectrum thus produced, was shown, in the 
energy ranges concerned here,* to remain con- 
stant at accelerating voltages equal to or greater 
than the energy of this isochromat. Hence, one 
should expect the activity in a sample irradiated 
with thin target x-rays to rise quickly to a certain 
value when the threshold energy for the process 
is reached and then to remain essentially constant 

, h > W. Barnes and P. W. Aradine, Phys. Rev. 55, 50 
ad Ie Lark-Horovitz, J. R. Risser, and R. N. Smith, Phys. 


Rev. 55, 878 (1939). 
7G. B. Collins and B. Waldman, Phys. Rev. 57, 1088 


(1940). 
8 E. Guth, Phys. Rev. 59, 325 (1941). 
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Fic. 1. Thick x-ray t arrangement with silver counter 
in position. 


until a second activation level is reached, where- 
upon the activity will again rise in a similar 
manner. The positions of these sudden rises in the 
activity vs. accelerating potential curves then 
indicate the energy of the activated levels. 
Similarly, as has been pointed out by Guth,® 
if we consider the thick target continuous x-ray 
spectrum as composed of a multitude of thin 
targets superimposed one behind the other, the 
activity should rise linearly with accelerating 
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potential until another activation level is reached, 
whereupon the slope of the activity curve in- 
creases. The intersection of the two linear por- 
tions thus indicates the position of the activation 
level. 

If the process is one of resonance absorption, 
the electron excitation curves should show sharp 
peaks when an activation level is reached and 
then fall to a lower value once this energy has 
been exceeded, similarly to the well-known exci- 
tation curves of Franck and Hertz for electrons 
striking an atom. 


EXPERIMENTAL PROCEDURE 
A. Silver 


It had previously been found® that silver pos- 
sesses a metastable state 93 kev above the ground 
level which decays with a half-life of 40 seconds. 
It has also been shown that this activity can be 

*L. W. Alvarez, A. C. Helmholz, and E. Nelson, Phys. 
Rev. 57, 660 (1940); G. E. Valley and R. L. McCreary, 
Phys. Rev. 56, 863 (1939); A. C. Helmholz, Phys. Rev. 60, 


oe get) ; E. Segré and G. T. Seaborg, Phys. Rev. 59, 212 


x103 


coun 
2 min 


activity in 


1.1 1.5 


3.5 


2.3 2.7 3.1 


beam energy in mev 


Fic. 2. Thick target x-ray excitation curve for silver. The positions of sudden changes in slope give the position 


of the energy levels which combine with the metastable state. 
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produced by direct x-ray bombardment of silver.!° 

In the present work, the x-rays were produced 
by a beam of electrons impinging on a gold 
target. The electrons were accelerated by a high 
pressure van de Graaff-Herb generator."' The 
sample of silver being irradiated formed the 


MEV SPIN 

3.13 

2.76 

2.32 

1.95 

1.59 

1.18 5/2 

0.093 riviy 
0 


Fic. 3. Nuclear energy level diagram for silver showing 
the various transitions by which nuclei may reach the 
metastable state. 


cylinder of a counter which was placed directly in 
. front of the target as shown in Fig. 1. 

The counters used were self-quenching and 
were filled with an argon (90 percent)—ether (10 
percent) mixture. It has been found that such a 
counter could be irradiated for several hours with 
very intense x-ray beams without showing a 
change in threshold or counting efficiency. Argon- 
alcohol counters, on the other hand, quickly 
deteriorate under the action of an intense x-ray 
beam. 

In obtaining an excitation curve it is necessary 
to find the activity produced as a function of the 
accelerating voltage, keeping all other quantities 
which might influence the activity constant. For 


( 10]. Feldmeier and G. B. Collins, Phys. Rev. 59, 937 
1941). 

" This name seems appropriate for generators of the high 
pressure type since much of the early work in their develop- 
ment was done by Herb and his group. 
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the work with thick target x-rays, the electron 
beam used was 200 microamperes, the irradiation 
time was fixed at two minutes, and the activity 
was taken as the number of counts above back- 
ground during the two-minute period starting 
five seconds after the beam was turned off. 

The curve thus obtained is shown in Fig. 2. It 
can be seen that the threshold, which represents 
the first activation level above the metastable 
state, is found at 1.18 Mev and higher activation 
levels, indicated by the sharp breaks in the curve, 
are found at 1.59, 1.95, 2.32, 2.76, and 3.13 Mev, 
respectively. 

‘From these data it is possible to draw a nuclear 
energy level diagram for silver (Fig. 3) indicating 
the excitation and the subsequent decay of the 
excited levels combining with the metastable 
state. 

Difficulties were encountered in observing the 
excitation of silver by direct electron bombard- 
ment since the energy of the conversion electrons 
emitted in the transition from the metastable to 
the ground state is very low (66 kev). These 
electrons are completely stopped by a very small 
thickness of material. As a consequence, it was 
necessary to find an arrangement by which the 
silver could be bombarded and then placed 
directly inside of a counter with as short a time 
delay as possible. To accomplish this, a counter 
was mounted on a large brass stopcock (Fig. 4) 
along with a ballast bottle filled with an argon- 
ether mixture at the same pressure as that in the 


counter. With this arrangement, the sample was | 


placed in a cavity (a) in the inner part of the 
stopcock when at position (0). After the cavity 
was evacuated at position (c) by means of a fore- 
pump, it was turned to position (d). The cavity 
was then filled with the gas mixture to the same 


liter 
ballast 


Fic. 4. Geiger-Miiller counter mounted on a brass a 
cock for admitting active sample to the inside of the 
counter. 
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Fic. 5. Thick target x-ray excitation curve for cadmium. 


' 
pressure as the counter after which it was turned 
to position (e) within the counter. With this 
arrangement the sample could be gotten into the 


counter within two minutes after irradiation was 
stopped. Activity produced by the electron bom- 
bardment was observed but the complete excita- 
tion curve has not as yet been studied for this 
element. 

The half-life of the metastable state was de- 
termined from a decay curve extending through 


six half-life periods with an initial activity of 8000 
counts/15 seconds. These data give a value of 
40.4+0.2 seconds for the half-life. 

This metastable state may belong either to” 
or both of which have spins of in 
the ground state. From the decay period and 
energy of the metastable state, this would imply 
a spin of 9/2 for the metastable state” and a 
probable spin of 5/2 for the level at 1.18 Mev. 


2 A. C. Helmholz, Phys. Rev. 60, 415 (1941). 
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Fic. 6. Electron excitation curve for cadmium. The positions of sharp rises give the 
energy of the excited levels in the resonance absorption. 
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Fic. 7. Absorption curve for the conversion electrons 
emitted in the transition: Cd*-+Cd++y. 


B. Cadmium 


The fifty-minute period of cadmium’ has 
been studied in a manner similar to that used 
with silver. The positions of the activated levels 
have been found from a thick target x-ray exci- 
tation curve. The activity was taken as the 
number of counts obtained during the first half- 
hour after irradiating a counter with a cadmium 
cathode for a time of ten minutes with a beam 
current of 250 microamperes. This curve is 
shown in Fig. 5 and as is seen, the threshold for 
this reaction, Cd+y—Cd*, is 1.25+0.03 Mev 
and higher activation levels are found at 1.68, 
2.08, and 2.56 Mev, respectively. In the region of 
the threshold, the samples were irradiated for 
much longer times with beam currents up to 700 
microamperes so that the activity of the cad- 
mium would greatly exceed the normal back- 
ground of the counter. These activities were then 
corrected to the proper value of irradiation time 
and beam current. 


An electron excitation curve for cadmium was . 


also obtained by bombarding thin samples of 


138 M. Dode and B. Pontecorvo, mag rendus 207, 287 
ter M. L. Wiedenbeck, Bull. Am. 


hys. Soc. 19, No. 3 
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cadmium directly with the electron beam. The 
activity produced in the samples was measured 
by the counter previously described (Fig. 4). The 
excitation curve (Fig. 6) shows a series of rather 


sharp peaks superimposed on a continuous back- 
ground. The positions of these peaks occur at the 


same voltages as the breaks in the thick target 


x-ray excitation curves. The relative sharpness of 
these peaks in the electron excitation curve gives 
further evidence that the process is one of line 
absorption. 

The energy of the metastable state was also 
determined by measuring the absorption of the 
conversion electrons in aluminum foils. Several 
runs were taken at a constant accelerating voltage 
and with constant beam current with various 
thicknesses of aluminum between the cadmium 
and the counting region. The activity was thus 
found as a function of the thickness (Fig. 7). The 
energy of the metastable state is then the energy 
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Fic. 8. Nuclear energy level scheme for cadmium. 


of the conversion electrons plus the binding 
energy of the K electrons of cadmium. The value 
thus found is 195 kev. 
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From these data one can draw a nuclear energy 
level diagram for cadmium (Fig. 8) showing the 
various energy levels combining with the meta- 
stable state. 

The isotope of cadmium which possesses the 
metastable state can be either Cd!® 112,118 since 
this activity has also been produced by Scherrer 
et al.* by the Cd reaction. Of these, 
can probably be eliminated in as much as ma- 
terial with this period cannot be grown from 
silver having the 22-second period. 


DISCUSSION 


It is of interest to notice that metastable levels 
were found for two “triads:’’ silver, cadmium, 
and indium and platinum,'* gold,’? and mer- 
cury.'® It is very likely that the group, kryp- 


“OQ. Huber, O. Lienhard, P. Scherrer, H. Waffler, Helv. 
Phys. Acta 16, 228 (1943). 

16 B. Waldman and M. L. Wiedenbeck, Bull. Am. Phys. 
Soc. 17, No. 5 (Nov. 1942); M. L. Wiedenbeck, Thesis 
on nuclear excitation of Indium. Both thick and thin target 
x-ray excitation curves and the electron excitation curve 
has Leen obtained for indium and will be reported in the 
near future. 

1 R. Sherr, K. T. Bainbridge, and H. H. Anderson, Phys. 
Rev. 60, 473 (1941). 

17 Data on gold will be published in a forthcoming issue of 
Phys. Rev. 

18 F, A. Heyn, Nature 139, 842 (1937); E. McMillan, M. 
Kamen, and S. Ruben, Phys. Rev. 52, 375 (1937); M. L. 
a. M. Cork, and R. L. Thornton, Phys. Rev. 52, 239 
(1937). 


ton,'® strontium,?® and rubidium, may form 
another such “triad.” 

The level schemes in the “triad,” silver, 
cadmium, and indium, show that the spacing of 
these levels is approximately uniform in any one 
of the elements. A similar condition has been 
found to exist in the case of nitrogen.” However, 
it should be noted that nitrogen is among the 
light elements while the elements studied here are 
of intermediate weight. It is possible that with 
the extension of these excitation curves to higher 
energies the levels will converge as in the case of 
atomic spectra. 
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An approximate expression for the partition function in a generalization of Bethe’s theory of 
superlattices including long range interaction is obtained as an integral of the approximate | 
energy expression in the generalized Bethe’s theory. An alternative form of the energy expres- 
sion is also considered, and higher approximations are treated. The paper ends with an applica- 


tion to the problem of adsorption. 


1. INTRODUCTION 


~ CE the appearance of Bragg and Williams’ 

first paper many methods for the approxi- 
mate solution of the problem of superlattices 
have been devised, of which the method of Bethe 
is quite a good one. In Bethe’s method the 
equilibrium values of the superlattice order and 
other quantities are obtained by an indirect 
method without evaluating the partition func- 
tion. In some problems, such as the critical 
condensation phenomenon in adsorption, the true 
critical temperature in a face-centered cubic 
lattice, and the separation into phases studied by 

_ Easthope,? the value of the partition function 
itself is indispensable for their investigation, and 
Bethe’s theory becomes powerless here. 

Some time ago the author put forward a theory* 
of obtaining an approximate expression for the 
partition function in Bethe’s theory, in its 
generalized form including long range inter- 
actions. Upon closer investigation it is found that 
the theory has two defects. First, the theory as- 
sumes that Bethe’s method is consistent when 
applied to the case of neighbor interaction so that 
the condition of integrability is satisfied and a 
partition function always exists. Although this is 
true in the cases considered then, it is found that 
in the case of Bethe’s second approximation and 
in the case of face-centered cubic lattice this is 
not true. Secondly, the numerical work involved 
in the calculation of the partition function is 


* At present a member of National Southwest Associated 
University, Kunming, China. ; 

1 See sg en review by F. C. Nix and W. Shockley, 
Rev. Mod. ae 10, 1 (1938). 
*C. E. Eas ope, Proc. Camb. Phil. Soc. 34, 68 (1938). 

+ J. S. Wang, “Free energy in the statistical theory of 
order-disorder transformation,’ Science Report of National 
Tsing Hua University [A], IV, Nos. 4-6 (April, 1941), 
(printed but failed to appear). 
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too laborious to be of any value in practical 
application. 

In this paper a modification of the theory is 
given. The present theory is perfectly general and 
applicable in all cases. Moreover, the numerical 
work has been lessened considerably so that 
practical applications are possible. The general 
theory is given in Section 3 and Section 4 below. 

In the application of Bethe’s theory one works 
with what may called the local grand partition 
function, instead of the ordinary partition func- 
tion. The use of the grand partition function in 
the theory of superlattices has been justified in a 
general way in the paper of the author cited 
above. It is found, however, that a slight modifi- 
cation of the theory there given helps to clarify 
some points concerned with the procedure of 
obtaining the approximate partition function 
from Bethe’s theory. This is done in Section 2. 

In the problem of adsorption the number of 
adsorbed atoms depends on the presence of the 
gas from which the adsorbed atoms come, and it 
is evident that the grand partition function 
should be used. But since the grand partition 
function is also used in the theory of superlattices, 
the theory of adsorption is mathematically the 


same as the theory of superlattices. The corre-. 


sponding quantities in the two theories have, 
however, different meanings, and some of them 
have to be treated differently. Thus the problem 
of adsorption receives a special treatment in 
Section 5. 

The notation used in the present paper is 
nearly the same as that employed in the previous 
one. To facilitate the reading of the paper it is 
perhaps better to explain some of the basic 
quantities first. We consider a binary alloy com- 
posed of two kinds of atoms A and B witha total 
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number N. The number of A atoms will be de- 
noted by N@, and that of B atoms by N(1—8). 
The N sites occupied by the N atoms are divided 
into two classes a and 8, the number of a-sites is 
Nr and that of 8-sites N(1—r). The value of r is 4 
for a simple cubic and a body-centered cubic 
lattice, and is } for a face-centered cubic lattice 
(r=} is also possible, e.g., AuCu). To simplify 
the writing we sometimes use the abbreviation 
N, for Nr, and Ng for N(i—r). In any given 
state of order let there be Nré, A atoms in a-sites 
and N(1—r)6s A atoms in £-sites. Since the total 
number of A atoms is N@, we have 


(1) 
The degree of order s is defined by 
(2) 


This definition of s reduces to that of Bragg and 
Williams in the special case of 0=r. 

Other symbols will be explained when they are 
first introduced. 

Recently a very interesting generalization of 
the quasi-chemical method of Fowler and 
Guggenheim has been made by Mr. C. N. Yang.‘ 
This method is very powerful, and perhaps 
Bethe’s method will be superseded by it. The 
chief merit of the method is that the approxi- 
mate partition function can be easily obtained as 
a closed expression. Mr. Yang is now making 
further investigation of the method. 


2. THE GRAND PARTITION FUNCTION 


The (configuration) energy of the alloy in the 
general case of long range interaction has been 
found to be® 


E= Z,( V,45— V,*) (3) 
= > M,44 V,, V, = V,44 + 2 


where M,44 is the number of AA pairs of atoms 
situated at a distance a, apart, v=1, 2, 3, ---; 
V,44, V,28, V,4® are the interaction energies of 
AA, BB, AB pairs at the distance a,; Z, is the 
sum of z, for all sites of the crystal, z, being the 
number of sites situated at a distance a, from a 


‘C. N. Yang, Chinese J. Phys. (in press). . 

‘Equation (3) is given in J. S. Wang, reference 3, 
Eq. (36). It can be easily obtained by the procedure given 
in Ng Wang, Proc. Roy. Soc. A168, 56 (1938), especially 
p. 95. 


given site. The sign > in (3) denotes a summa- 
tion over the values of ». 

From the above expression for the energy the 
ordinary partition function is obtained as 


Po=exp 
X {exp [— — 
XLOQ(Ga, %,T), (4) 
Q(6a, T) = exp (—M,44V,/kT), 


where ym is the number of arrangements for given 
values of 6., 63, M,44. The summation of 
Q(6., 93, T) in the expression for P, is taken over 
all 6, and 6g satisfying the relation (1) with a 
fixed @, and the summation in the expression for 
Q(6., 93, T) is taken over all possible values of 
M,44 for given and 

The purpose of the present paper is to obtain 
an approximate expression for Q(6., 63,7) in 
Bethe’s method. 

The equilibrium values of 6., 63, and E; are 
given in terms of Q by the usual formulae 


or 
a 
E,=kT*— log 0. 6 
1 (6) 


These equations can be deduced from the stand- 
ard formula in terms of Ps on replacing P, by its 
maximum term Q. 

Let us now form the grand partition function 
K by considering the number of A atoms as 
variable. We have 


K=>\"*P,=exp V,28/kT) Ki, 
Ki =D (0., Os, T), (7) 
t=) exp [— DZ,(V,48— V,28)/NkT], 


where the summation in K, is taken over all 
possible values of @, and 6s. In order to obtain the 
equilibrium values of @., %, E; from the grand 
partition function, it is more convenient to 
introduce two ¢’s, & and és, say, such that K, 
becomes 


Ki= 6s, T). (8) 


= 
} 
| 
) 
| 
| 
is 
ind 
cal 
hat 
Ow. 
rks 
ion 
nc- 
in 
na 
ted 
lifi- 
rify 
of 
‘ion | 
of 
the 
d it 
tion 
tion 
ces, 
the 
ive, 
1em 
lem 
in 
r is | 
ious | 
it is 
asic 
om- 
otal | 


Then the usual formulae for the equilibrium 
values of 0., 0s, E; are 


Nré,= log Ki, 

(9) 
N(i—r)6s= log Ki, 


E,=kT*— log 10 
1 aT og Ai (10) 


and in the final result we must put 
fa= (11) 


To prove that the equilibrium values deter- 
mined from the grand partition function K, by 
means of (9)—(11) are the same as those calcu- 
lated from the ordinary partition function Q by 
means of (5) and (6), we may adopt the procedure 
of replacing K, by its maximum term. The 
maximum term in the summand of K, given by 
(8) occurs at the values of 0, and 6 satisfying 


1 @ 
og & = 0 
(12) 


log No 305 Q=0. 
When K; is replaced by its maximum term, 
Eqs. (9) become identities, Eq. (10) reduces to 
(6), and Eq. (11) reduces to (5) with the help 
of (12). 

Thus the use of grand partition function is 
justified. 

A more rigorous proof of the above result by 
means of contour integrals has been given in the 
paper of the author cited above. This proof will 
be briefly indicated here with a slight change. 
Consider the series 


S(x,y) = Os, T). 


Thus K, given by (8) is S(é., &s) and that given 
by (7) is S(é, —). By Cauchy’s theorem we have 


S(x, y) 


93, T)= 


(2m)? at 


By the general theory of statistical mechanics 
given in Fowler’s book, the integrand has a 
unique maximum modulus at the saddle point 
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x=£,, determined by 


Na log S(Eas fs), 


0 
log S(&a, &). 


This proves (9). Differentiating Q with respect to 
6. we have 


IO 


There is a similar equation for 0Q/063. Using the 
method of steepest descent we obtain (12). 

To determine the equilibrium value of 6, for a 
given value of @, we put y=x in S(x, y) and obtain 


S(x, 


S(x,y) 


—dx 


and 
NaGaQ(Ga, 98, T) 
S(x, x) 


1 
=— — log S(x, ———dx, 


where the symbol >-’ denotes summation over 6, 
with 63 given as a function of 6. by (1). Hence the 
equilibrium value of N.0, is given by 


where ¢ is the saddle point determined by 


= log S(é, ). 
N@ (é, &) 
This proves (11). 


3. APPROXIMATE PARTITION FUNCTION 


In Bethe’s method one fixes one’s attention on 
a group of quite a small number of sites and 
constructs the grand partition function for these 
sites; the effect of omitting other sites of the 
crystal is taken account of by introducing suitable 
parameters for the outermost sites of the group. 
The group of sites are usually divided into the 
central site, the first shell sites which are nearest 
neighbors to the central site, the second shell 
sites which are nearest neighbors to the first shell 
sites, and so on. In the first approximation the 
group of sites considered include only the central 
site and the first shell. 
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The extension of Bethe’s first approximation to 
include long range interaction has been given in 
an earlier paper,® and the results will be quoted 
- here. The equilibrium values of 0, and 6% are 
given by 


8a 


exp (— TI, 


(13) 


68 G, 
—U5/kT)—, 
3 exp (— 


where F and G are functions of the parameters 
éx, €s, and £4, respectively, and the subscript 7 
attached to a function means that the variables 
éz, etc. in that function are changed to mea, etc., 


where 
=exp (— V,/kT), 


The parameters €., €3, fa, (3 are determined by the 
following equations: 


= (1— Oa) 
Zapp = (1 — Oa) 
Zpa9a= (1— Ogu», 
= (1 —03)x 


(14) 


(15) 
where 
e=ta—log F, p=es— log F; 


0 
= ‘a— | G, = ] G. 
w og x og 


The quantities U, and Us are given by 
V1) (Uap — V1) 0p, 
(Ups — Vi) Op. 

The equilibrium value of the energy £; is given 
by 

Vit Ue} 

(17) 


In order to obtain an approximate expression 
for the partition function Q consistent with 
Bethe’s approximation, it seems that a natural 
way would be to assume that the equilibrium 
values of 62, 93, and E,; given by (13) and (17) are 
correct, and to substitute the values of & and & 
given by (13) into (12) and the value of E; given 
by (17) into (6) to obtain Q by integration. 


(16) 


§J.S. Wang, Proc. Roy. Soc. A168, 56 (1938). 
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Unfortunately this procedure does not work, be- 
cause the system of differential equations for Q 
so obtained does not, in general, satisfy the con- 
dition of integrability. 

In my previous paper this difficulty was avoided 
by regarding U. and Us as unknown functions to 
be determined by the condition of integrability. 
Asystem of partial differential equations for these 
two quantities was thus obtained, whose solu- 
tions were made unique by the condition that U, 
and Us reduce to the expressions given by (16) in 
the special case 0.=63 and in the special case 
@=1. These conditions are evidently obtained 
from physical requirement. 

The physical basis for the above procedure is 
that the assumption of a uniform distribution of 
A atoms with a fraction @, on a-sites and a 
fraction 63 on 8-sites (outside the small group of 
sites) used in deriving (16) is not true and should 
be replaced by the condition of integrability. But 
the assumption is certainly true in the special 
cases 0,=63 and @=1, and so (16) serves as a 
boundary condition. 

It has been mentioned in the Introduction that 
the above theory suffers from two defects. The 
most serious defect is that in the special case of 
neighbor interaction the two quantities U, and 
Us should be zero, but then the condition of 
integrability is not satisfied in general. If we 
again introduce U, and Us into the formula to 
restore the condition of integrability, it is hard to 
interpret their physical significance. A second 
defect is that the numerical work involved in 
calculating U, and Us is too laborious. 

We shall now modify the theory by preserving 
formulae (16) and (17) for U., Us, and E,, but 
regarding £ and é given by (13) as incorrect. We 
shall obtain Q by integrating (6), and take (12) as 
the equations determining & and ts. To integrate 
(6) we remember that as T— ~ , Q approaches the 


value 
83) = ( ( (18) 


Hence 


log Q=1 (19) 


in which the value of E; given by (17) should be 
used. 
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When Q has been obtained in this way, the -site. Thus 
equilibrium values of @, and 6% will be given by = ” 
(5), instead of by (13) and (11). (20) 

We shall now transform the integral in (19) to The quantities Aa,» €a, etc., are ‘determined as 
a simpler form. Let us denote the local grand functions of 02, 5, T by (13)-(15). As T+, 
_partition function for central a-site by K., and these quantities reduce to 
write A, in place of £ exp (— U./kT), remember- a 6s 
ing that has no real connection with £, in (12) » (21) 
owing to the breakdown of the integrability con- , ° 
dition. Similarly, we have Kg and g for central Then the result of integrating (19) is 


log Q=log Ka— 0a log €a — log €3} 
{log 0s log Xs — log fa — 2990s log 
log log (1 log 9) log (1 — 4s) 


Vv; 
+4f { Na(1 — 0a) F*+ (22) 


where F* is defined by 
_kT’ log Flee, (23) 


Since it follows from Stirling’s formula that 
log g= — Na[0a log log (1—0a) ]— Op log O-+ (1 — log (1— 4) (24) 


we have 
log ]+4N log G,+(1 — Og) log G— log fa — log 


d 
+} f { Na(1 — 02) F*+ Nee F,*+Na(1 (25) 


The last integral should be evaluated by numerical integration. 
4. HIGHER APPROXIMATIONS 


The method of obtaining the partition function in the first approximation given in the preceding 
section is evidently applicable to higher approximations. The group of sites considered in the general 
case is divided into the central site, the first shell sites which are at a distance a, from the central site, 
the sites which are at a distance a2 from the central site, and so on. Let the numbers of a- and -sites 
among those which are at a distance a, from the central site be 2,aa, 2rag if the central site is an a, and 
” Zy8ay 26g if the central site is a 8. These numbers satisfy the relations 


Let us consider the group with a central a-site. We shall introduce a parameter dq for an A atom on 
the central site, and parameters ¢,. and ¢,s, respectively, for an A atom on the a- and §-sites which are 
at a distance a, from the central site. In any distribution of A atoms in the group of sites considered, 
let there be m) A atoms in the central site (n9=1 or 0 according as the central site is occupied or not), 
n,=n,' +n," A atoms in the sites at a distance a, from the central site, ,’ being the number in the 
a-sites and n,”’ the number in the §-sites of them. Let m, be the number of pairs of A atoms in the 


(26) 
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group at a distance a, apart, and let 7. be the number of arrangements for a given set of values ,’, 
n,'’, and m,. Then the local grand partition function is _ 


Ka=F+aFy 
with ” 


F, = 1€18, N2€28, 
in which the suffix « denotes the outermost sites of the group and the relations 


m,=non,+m,’, (vSx) 
have been used. 
The parameters will be determined by the conditions 


0 
0a=a— log Ka, 


These equations will also be expressed in a new notation, analogous to (14), 


= (1 0a) Ba Pras 98 (1 0a)Wr+ 
where 


log F, log F. 
0,8 

There are a similar set of equations for the group with central £-site. 

The present formulation of the higher approximations is different from Bethe’s original formulation 
and from the modified one given by the author’ for the case of adsorption with long range interaction. 
Bethe considered only the special case of equal A and B atoms and treated them symmetrically, so 
that in the first approximation he put A, =1 and in the second approximation he put A2=1, ¢s=1 (in 
Bethe’s case there are no a-sites in the first shell surrounding a central a-site). When the numbers of 
A and B atoms are not equal, it is more convenient to concentrate attention solely on A atoms and 
to introduce a A, different from unity, as first done by Peierls.* Such a procedure also works in the case . 
of equal A and B atoms. Peierls put \.= a= £ as the condition of equilibrium. This was also followed 
by the author®® in the case of long range interaction by putting \.=¢ exp (— U./kT), the second 
factor being introduced to take account of the influence of exterior sites. Since this does not lead to 
consistent results when the partition function is sought, we have replaced it by Eq. (5) in the preceding 
section. In accordance with a procedure adopted in the theory of adsorption with long range inter- 
action,’ a natural generalization to higher approximations would be to put 


exp (—Ura/kT), exp (30) 


where U;, is the energy owing to exterior sites on an a-site at a distance a, from a central a-site, 


and U,. is the energy owing to exterior sites on a §-site at a distance a, from a central a-site. Such 
a generalization evidently fails for the same reason as Peierls’ assumption on A, does. Therefore 
we hana it by the condition (28,). 


LS W Wang, Proc. Camb. Phil. Soc. 34, 238 (1938). 
Peierls, Proc. Roy. Soc. A154, 207 (1936). 
*See J. S. Wang, Proc. Roy. Soc. A168, 68 (1938). 
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The approximate partition function will still be given by (19), in which the energy now has the 
form 


where U, and Us are given by 


(32) 
If we introduce the notation 
F( (= VAL v.) Tes ) II (33) 
V; aT og Ela, €18, E2a, Ya m, m, €va 


similar to (23), then 


log Q=log e-1( 145 log 


foe x. — a log Xa Zraa9a >, 


1/T Vi 1 


The actual calculation may sometimes be simplified by assuming reasonable values for some of the 
a, etc., and dropping the corresponding equations in (28,). Thus in the case of the second approxi- 
mation, for example, one may assume €i¢ and ¢€:g to be given by (30), with £4, & determined by the 


first approximation, and drop the corresponding equation (28,). In a certain sense, it is better to drop | 


the equation (28,) and retain (28;) to determine €2¢ and €2. This alternative procedure has the 
advantage in making the approximation better for the central site and the first shell, and actually 
corresponds to Bethe’s original procedure. Bethe also assumed one and the same parameter for A 
atoms in the corner and medium sites of the second shell. This procedure has the disadvantage of 
losing the symmetry, and it is perhaps still better to retain all the equations (28,), in which the corner 
and medium sites are distinguished by the difference in », i.e., in their distances from the central site. 

A different way of calculating the energy has been adopted by Mr. Yang in his generalization of 
the quasi-chemical method. His method can also be slightly modified to apply to the present problem. 
In obtaining the expression (31) for the energy we have assumed that the equilibrium values of M,44 


for yx are given by 
(mony) 


We shall now modify this assumption for y= 1 but shall retain the assumption for vy>1. In the case of 
v=1, let us assume that m, is the same fraction of u as (M144), is of $Z2,;=4Nz, where u and 4Nz are 
the total numbers of pairs of neighboring sites for the group and for the whole crystal, respectively.” 


1° The value of u is the same for central a- and for central 8-site on account of equal z for them. 
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Under this assumption, we have 


oT K. 


Summing over all central a-groups and all central 6-groups we obtain 


Nek log Ka log Ks= Ne (X )atNe(X m,V,)s. 


Hence the modified energy expression is 


y>1 vy>1 


»>1 


When this expression is substituted into (19), the partition function becomes, after integration, 


log Q=log log 


+—N. l | 
{log 


+7 Nollog log log fra— 


— 0a log lOg log €,g} 


LX log 


2uk »>1 
is V, 


WT 


1 
{ N, Pry +N, 298 (wry + } d (=) (36) 


The advantage of this alternative formulation in the case of neighbor interaction is evident. For 
then the integral in (36) is identically zero and the partition function is obtained as a closed expres- 
sion. In the particular case of first approximation for the simple cubic and body-centered cubic 


lattices these two formulations are the same. 


It may be argued that since Bethe’s method would give a better equilibrium value for the central 
‘site than for the outer sites, the expression (31) for the energy is perhaps better than (35). 


5. THE PROBLEM OF ADSORPTION 


In the problem of adsorption with a single kind 
of adsorbed atoms, the theory of superlattices can 
be applied provided the A atoms are interpreted 
as the adsorbed atoms. It is now more convenient 
to choose 6, s, T as independent variables, the 
relations between 6, s and 6., 63 are given by (1) 
and (2). The equilibrium value of s as a function 
of @ and T is still given by (5). But @ depends on 


the pressure p of the gas, and so instead of 
regarding (12) as equations determining & and 
t; in terms of @ and s, we regard them as de- 
termining in terms of and T, with connected 
with the pressure by the relation :" 


hbex!*T v.(T) 37 


1 See J. S. Wang, Proc. Roy. Soc. A161, 127 (1937). 
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In (12) becomes 


1 
log tog Q. (38) 
This equation gives the relation between @, p, 
and 7, and is commonly known as the adsorption 
isotherm. 

The heat of adsorption g is defined as the 
decrease of energy when one atom is transferred 
from the gas phase to the adsorbed phase and is 
given by" 

1 dE; 


N 00" (39) 


d 
qo=x+ kT +kT?— log —, 
dT 


and the equilibrium value of s determined by (5) 
should be substituted in it. 

The special case of equal numbers of a- and 
_ B-sites is interesting for its simplicity. It is easy 
to see that in this case EZ; is symmetrical in a 
and 8, and so it is an even function of s. To study 
the nature of the dependence of E; on 6, we shall 


consider the first approximation for simplicity. It ~ 


will be evident that the higher approximations 
can be treated in the same way. In the case of 
first approximation, there is a single « and a 
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where £; is considered as a function of @ and T 


single ¢. The functions y and w satisfy the 
relation® 

(f=yorw). (40) 
It then follows from (14) and (15) that the 
transformation 


(41) 
changes 0, to 1—6g and @ to 1—@q, so that s is 


unaltered while 6 is changed to 1—8@. Then (17) 
shows that 


is an even function of @—4. Consequently 


U 
log and g—gqo+}U 
are both odd functions of @— 4 and even functions 


of s. 

The same conclusion holds in the case of 
vanishing order (s=0) in a general lattice (i.e., r 
need not be 4). For then 0,.=6s=0, ea=es=f, 
={g=e, Fe, =G(e, €-)=f(e), and in the first 
approximation becomes ; 

(42) 
It is easily seen that E,—4NU(@—}4) is an even 


function of @—4, and so the same conclusion on ¢ 
and g follows. 
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The Multiple Scattering of Waves* 


I. General Theory of Isotropic Scattering by Randomly Distributed Scatterers 
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While the problem of the multiple scattering of particles by a random distribution of scat- 
terers has been treated classically through the use of the Boltzmann integro-differential 


equation, the corresponding problem of the multiple scattering of waves seems to have received 
scant attention. All previous treatments have considered the problem in the “geometrical 
optics” limit, where the rays are regarded as trajectories of particles and the treatment for 
particles is then applied, so that the interference phenomena in wave scattering are neglected. 
In this paper the problem of the multiple scattering of scalar waves by a random distribution 
of isotropic scatterers is considered in detail on the basis of a consistent wave treatment. The 
introduction of the concept of “randomness” requires averages to be taken over a statistical 
ensemble of scatterer configurations. Equations are derived for the average value of the wave 
function, the average value of the square of its absolute value, and the average flux carried by 
the wave. The second of these quantities satisfies an integral equation which has some simi- 
larities to the corresponding equation for particle scattering. The physical interpretation of 


the results is discussed in some detail and possible generalizations of the theory are outlined. 


INTRODUCTION 


HE problem of the multiple scattering of 


particles by a random distribution of scat- 
terers has been treated in some detail in litera- 
ture.! These treatments are all based essentially 
on the Boltzmann integro-differential equation 
of conservation in phase space and are thus 
classical rather than quantum-mechanical in 
foundation. On the other hand, a consistent 
treatment, based on the wave equation itself, of 
the multiple scattering of waves by a randomly 
distributed collection of scatterers, has not, to 
the writer’s knowledge, yet appeared in litera- 
ture. Those treatments which have appeared, for 
example, in discussions of the problem of the 
propagation of light through stellar atmospheres 
and through turbid media,’ are all based on the 


* Publication assisted by Ernest Kempton Adams Fund 
for E Physical Research of Columbia University. 
dence should be addressed to 64th Floor, 
350 Fifth Avenue, New York 1, New York. 
1E. Fermi, Ricerca Scienta VII-II, 13 (1936); Yost and 
Dickinson, Phys. Rev. 50, 128 (1936); G. C. Wick, Atti 
del Acad. Reale dei Lincei 23, 775 (1936); L. S. Ornstein, 
Kon. Akad. van Wet. te Amst. Proc. XXXIX, No. 9 
(1936); Halpern, Lueneburg, and Clark, Phys. Rev. 53, 173 
(1938) ; Schuchard and Uehling, Phys. Rev. 58, 611 (1940). 
°E, Hop f, “Problems of radiative equilibrium,” Cam- 
bridge Tract No. 31 (1934); E. A. Milne, “Thermodynamics 
of stars,” Handbuch der Astrophys., Vol. III/I (Julius 
, Berlin, 1930); G. Jaffe, Ann. 195 
990); S. Q. Duntley, J. Opt. Soc. Am. 33, 252 (1934); 
Hulbert, J. Opt. Soc. Am. 33, 42 (194 
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approximation of considering light rays as the 
trajectories of particles (photons), which are then 
treated on the basis of the Boltzmann equation. 
This, in a sense, corresponds to treating light 
scattering in the ‘‘geometrical optics limit’’ where 
interference phenomena characteristic of waves 
is neglected. Such a treatment might be expected, 
a priori, to be valid if the wave-length of the 
waves is small compared to the mean free path 
of the associated particles in the scattering 
medium, but the only satisfactory justification 
for such a treatment lies in showing that the 
wave treatment leads to the same results. 

The investigation of the multiple scattering 
of waves is not only important for the above 
reason but also for several others. The multiple 
scattering of sound waves by the water droplets 
of a fog or by other collections of small obstacles 
often occurs in cases where the interference 
effects may not be legitimately neglected, and 
thus this problem cannot be satisfactorily dis- 
cussed except by a wave treatment. Also in the 
scattering of electrons and other fundamental 
particles by aggregates of nuclei, atoms, or 
molecules in gases, liquids, and amorphous solids, 
the use of the classical particle mechanics may 
be questionable, and the quantum-mechanical 
treatment of this problem necessitates the inves- 
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tigation of the multiple scattering of waves. 
Finally, the connection between the so-called 
coherent scattering (leading to the macroscopic 
phenomena of refraction and specular reflection) 


_ and the incoherent scattering can only be appre- 


ciated through a consistent wave treatment of 
scattering from beginning to end. 

In this paper we shall develop a wave theory 
of multiple scattering for the case of scalar waves 
which are isotropically scattered by a random 
distribution of scatterers, the ‘‘randomness” 
being defined explicitly in the next section. We 
deal with abstract waves of unspecified type so 
that the treatment is applicable to any type of 
wave satisfying the scalar wave equation. The 
concept of randomness being essentially statis- 


tical, it is necessary to concentrate one’s interest 


in the average values of physical quantities, the 
averages being taken over a statistical ensemble 
of collections of scatterers. The quantities of 
greatest interest are the wave function itself, the 
square of its absolute magnitude, or its mean 
square value, and the flux of some quantity 
(energy, probability, etc.) carried by the wave. 

The determination of the average value of the 
wave function is practically identical with the 
already treated theory of the index of refraction 
of material media (other than crystalline media) 
for light,’ but is included in the present treat- 
ment since it forms a necessary introduction to 
the second problem of determining the average 
value of the square of the absolute value of the 
wave function. It is shown that the average 
value of the wave function satisfies the wave 
equation in a continuous medium in which no 
scatterers are present and in which the velocity 
of propagation is different from the velocity in 
the original medium in the absence of the scat- 
terers and is, in general, a function of position 
and complex rather than real in value. In other 
words, the incident wave and the scattered 
waves combine on the average to form a wave 
which travels uniformly without scattering at a 
different velocity from the incident wave and 


with attenuation. This wave shows refraction | 


and reflection phenomena and thus describes the 
coherent aspects of scattering, such as the 
specular reflection taking place at a boundary 


use Born, Optik (Julius Springer, Berlin, 1933), pp. 
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between a region where scatterers are present 
(in sufficiently great density) and a region where 
no scatterers are present. An example would be 
the reflection of light at a glass-vacuum interface, 
the atoms of the glass being the scatterers. 
The principal contribution of this paper, how- 
ever, lies in the treatment of the average of the 
square of the wave function, which is in general 
different from the square of the average value 
of the wave function. The difference represents 
what is usually referred to as the incoherent 


scattering, the term incoherent referring not to 


a lack of definite phase relationship between the 
incident wave and the wave scattered by a par- 
ticular scatterer, but to the statistical super- 
position of these scattered waves when propa- 
gated to the point of observation because of the 
“random”’ distribution of the scatterers. The 
average of the square of the wave function will 
be shown to satisfy an integral equation which 
has some similarities to the form of the Boltz- 
mann integro-differential equation governing the 
multiple scattering of particles. This integral 
equation is just as fundamental to the problem 
of the multiple scattering of waves as the Boltz- 
mann integro-differential equation is to the 
problem of the multiple scattering of particles. 
The physical interpretation of the integral equa- 
tion may be fairly readily recognized: The equa- 
tion states that the average of the square of the 
wave function at any point is equal to the square 
of the average wave function at the point plus 
contributions representing scattering from all 
other portions of the medium; these contribu- 
tions are proportional to the average of the 
square of the waye function at each point and to 
a scattering cross section per unit volume at 
each of these points. The value of the scattering 
cross section per unit volume at each point and 
the mode of propagation from each point to the 
point of interest are exactly defined in terms of 
the scattering properties of the individual scat- 
terers and the properties of the statistical 
ensemble of scatterers studied. In fact the 
integral equation could be derived from its 
a posteriori physical interpretation, but in that 
case the value of the scattering cross section 
per unit volume and the mode of propagation 
of the scattered waves would have to be intro- 
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duced ad hoc, so that such a derivation could 
hardly be considered rigorous or fundamental. 


STATISTICAL CONSIDERATIONS 


Since there is no inherent limitation in the 
theory restricting it to the case where all scat- 
terers have identical scattering properties, we 
shall assume that the scattering properties of the 
scatterers are determined by a scattering param- 
eter denoted by s. In some applications s may be 
a continuously varying physical parameter such 
as the radius of water droplets in the problem 
of the scattering of sound by a fog, or it may be 
simply an index identifying particular types of 
scatterers such as various types of nuclei present 
in the problem of the scattering of neutrons by 
an amorphous solid or by a gas. In what follows 
we shall assume it to be a continuous variable, 
but the transformation of the equations in case 
it takes on only discrete values involves only an 
obvious replacement of integrals by summations. 

If we have a collection of N scatterers and are 
given for each its position and the value of s 
which characterizes its scattering properties we 
shall say that we have a particular configuration 
of the scatterers. Now, in the problem of mul- 
tiple scattering of waves by a collection of 
amorphously distributed scatterers, one rarely 
has sufficient information to establish the micro- 
scopic configuration of the collection and one is 
therefore not so much interested in the values 
of physical quantities for a particular configura- 
tion of the scatterers as in the average value of 
these quantities, the average being taken over 
the ensemble of possible configurations of the 
scatterers consistent with what information one 
has as to the macroscopic state of the collection. 
Thus if we specify the positions of the scatterers 
by their position vectors fr, f2, ---fy, and their 
scattering properties by the values of the scat- 
tering parameters 5), S2, ---sy, a particular con- 
figuration of the scatterers will be given when 
all the r,’s and s,’s are’specified. We may then 
indicate the ensemble of configurations in which 
we are interested by a probability distribution 
function P(ri, fe, 51, Se, so that 


P(r, To, Tw, $1, 


(1) 
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represents the probability of finding the scat- 
terers in a configuration in which the first 
scatterer lies in the element of volume dr; about 
the point r,; and has a scattering parameter 
lying between s; and s,;+ds;, the second scat- 
terer lies in the element of volume dr, about the 
point rz and has a scattering parameter lying 
between sz and s2+dse, etc. To find the average 
of a physical quantity over the ensemble of con- 
figurations we then multiply its value for the 
configuration (ri, Te, ***Ew, 51, * Sw) by (1) 
and integrate over all values of the r;’s and s;'s 
accessible to the scatterers. 

In this paper we shall be interested only in 
probability distribution functions in which the 
probability that a particular scatterer is located 
in some volume element and has a value of s 
located within some range ds is independent of 
the locations and scattering parameters of the 
other scatterers and is the same for all scatterers. 
In this case we may write for the probability dis- 
tribution function 


P=(1/N)¥n(ri, Sw). (2) 


The average number of scatterers per unit 
volume in the neighborhood of the point r having 
scattering parameters lying between s and s+ds 
is then 


n(r, s)ds, 
and the total number of scatterers per unit 
volume in the neighborhood of the point r is 
= f n(x, s)ds, 


the integration being taken over all values of s. 


.Finally the normalization of the probability 


requires 
s)dsdr = N 


where V represents the total volume accessible 
to the scatterers. 

We shall call the average of a physical quantity 
over the ensemble of configurations a configura- 
tional average and shall denote it by enclosing the 
quantity which is averaged in angular brackets. 
Thus if f is a function of the r,’s and s,'s, its 
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configurational average is 


fo Sf fof 


P(r;, (3) 


If we wish to omit the integration over the 
positions and scattering powers of one or more 
of the scatterers we shall indicate this by sub- 
scribing the indices of these scatterers to the 
angular brackets; thus (f); indicates that the 
integration over r; and s; is to be omitted in (3). 


STATEMENT OF THE PROBLEM 


In this section we shall make a precise state- 
ment of the problem which we shall treat in this 
paper. Interest here will be restricted to the case 
of steady-state scattering of waves of a single 
frequency w so that the value of the scalar wave 
function at the point r at time ¢ can be repre- 
sented as y(r)e**‘. In the absence of scatterers 
from the medium, y¥(r) will then satisfy the wave 
equation 


where ko =w/co and is the wave velocity in the 
scatterer-free medium. Furthermore we shall 
consider here only isotropic point scatterers, 
that is, those which scatter spherically sym- 
metrical waves (s-wave scattering only). In this 
case the wave function in the presence of the 
scatterers will have simple poles at the positions 
of the scatterers and in the neighborhood of the 
jth scatterer will behave like 


A;exp [—tko|r—r;| ] 


|r—r;| 


We may then define the external field acting on. 


the jth scatterer as 


A; —tko|t—T; 


|r—r;| 


this field having no singularity at r=r;. We then 


characterize the scattering properties of the- 


scatterers by the relationship 
A (6) 
making the strength of the scattered wave from 
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a scatterer proportional to the external field 
acting on it. In the above g(s;, w) will be referred 
to as the scattering coefficient for the jth scatterer 
and will be abbreviated to g;. 

The problem which we shall consider in this 
paper is then the following: Given the function 
g(s,w) and the distribution function n(r, s) for 
the scatterers, and given also the form of the 
wave function yo(r) which is present in the 
medium in the absence of the scatterers (the 
incident wave), to find the configurational 
average of the wave function, y¥(r), and of its 
mean-squared value, |¥(r)|*, in the presence of 
the scatterers. 


FUNDAMENTAL EQUATIONS 


Consider a particular configuration of the 
scatterers. We then attempt to represent the 
value of the wave function at the point r in the 
form 


+2 A;E(r, 
7 
where we have introduced the abbreviation 
exp [—iko|r—r’| ] 


E(r, 


The equation represents the field as the sum of 
the incident wave and spherical waves diverging 
from each of the scatterers. The external field 
acting on the jth scatterer is then 


= vol) + ty). 
We must then have, preneonas to Eq. (6), 
so our equations become 


(ry )E(r;, ty). (8) 


These represent the fundamental equations of 
multiple scattering. 

The direct method of solving our problem 
would then consist of solving the set of simul- 
taneous linear algebraic Eqs. (8) for the y*(r) 
and substituting these in Eq. (7), thus giving 
¥(r) as a function of the positions and scattering 
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parameters of the scatterers. Taking the con- 
figurational average of this quantity and of the 
square of its absolute value according to Eq. (3) 
would then give us the desired results. Unfor- 


tunately, it does not seem possible to carry this — 


procedure through because of the complexity of 
the necessary integrations, and it is necessary to 
resort to another procedure. This alternative 
method consists of attempting to find equations 
satisfied by (|¥(r)|) and (|¥(r)|*) and then 
solving these equations for the desired averaged 
quantities. We shall now carry through, with 
some approximations, this latter procedure. 


THE CONFIGURATIONAL AVERAGE OF THE 
WAVE FUNCTION 


Let us take the configurational average of 
both sides of Eq. (7) formally to obtain 


+E f f eXvi(rs)); 


XE(r, 


=Yo(r) + f G(r) (9) 
where 


G(r)= | g(s, w)n(r, s)ds. (10) 
The quantity (y*(r,;)); represents the external 
field acting on the jth scatterer averaged over 
all possible configurations of all of the other scat- 
terers. The only rigorous way of evaluating’ it 
seems to be to solve the set of simultaneous 
linear algebraic Eqs. (8) and then carry out the 
necessary integrations to obtain the above con- 
figurational average. This again is not feasible. 
We therefore must resort to the approximation 
of replacing the external field acting on the jth 
scatterer averaged over all configurations of the 
other scatterers by the average field which would 
exist at the position of the jth scatterer when the 
scatterer is not present. This last average field 
would differ only by a term of the order 1/N 
from‘ (y(r;)), so that if N is large, and if the 
above approximation is valid, which it appears 

‘¥(t;)=Lim y(r) is singular, but (¥(r;)) =Lim (y(r)) is 


tr; 


regular. 


f = 
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to be on physical grounds, we may substitute 
(v(r,;)) for (y4(r,;)); under the integral sign in the 
above equation thus obtaining the integral 
equation 


= volt) + f G(r’) E(x, 
Vv 


for the configurational average of the wave 
function. 

To see the physical significance of this equa- 
tion, let us apply the operator 


V+ke (12) 


to both sides of the integral equation, remem- 
bering that 
(V?+ ko?) E(r, r’) —4ri(r—r’), (13) 


where 6(r—r’) is the three-dimensional Dirac 
6-function defined by the equation 


f(r) if rlies in V 

0 if r lies outside of V. 

We then obtain 

= 
If we write this as 

=0, 
k*(r) =ko?+4nG(r), 


(15) 
(16) 


we see that (¥(r)) satisfies the wave equation in 
a ‘‘continuous medium” in which the velocity of 
propagation depends upon the scattering coef- 
ficients and density of distribution of the scat- 
terers and is, in general, a function of position. 
Thus the problem of finding the average value 
of the wave function has been essentially reduced 
to solving a boundary value problem in the wave 
equation. The boundary conditions are implied 
in the integral equation itself and depend on the 
function G(r). If G(r) is everywhere continuous 
and approaches a constant value or zero at 
infinity, then the boundary conditions are that 
(¥(r))—yYo(r) be everywhere continuous and 
have a continuous gradient and at infinity 
represent outward travelling waves. In another 
important case, that in which G(r) is sectionally 
continuous, the boundary conditions are that 
(¥(r)) —Yo(r) be everywhere continuous and have 
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a continuous normal derivative across a surface 
of discontinuity of G(r) and at infinity represent 
outward travelling waves, provided again that 
G(r) approaches a constant value or zero at 
infinity. In both cases, of course (y(r)) must 
approach y»(r) as G(r) becomes zero everywhere. 

This reduction of the determination of (y(r)) 
to the solution of a boundary value problem in 
the wave equation corresponds, of course, to the 
well-known theory of the refractive index for 
light of gases. ' 

We note also that the integral equation (11) 
can in principle be solved directly. One method 
of solution is the Liouville-Neumann method of 
successive substitutions or iteration method.® 
It consists in repeatedly substituting the ex- 
pression for (¥(r)) as given by the right-hand 
side of the integral equation for (y(r’)) under the 
integral sign, thus yielding in our case the 
infinite series 


= (17) 


where 


which, if it converges uniformly, is the desired 
solution. 

Revérting for a moment to the approximation 
((pi(r;))<~(WX{r,))) made at the beginning of this 
section, it should be mentioned that it would be 
desirable to demonstrate its validity mathemati- 
cally rather than on intuitional grounds. Unfor- 
tunately it has not yet been found possible to do 
this. A similar approximation must be made in 
the next section to which this remark also applies. 


THE CONFIGURATIONAL AVERAGE OF 
|? 


To obtain the configurational average of 
|¥(r)|? we follow an analogous procedure to 
that used for y(r), but the analysis now becomes 
somewhat more involved. To begin we first mul- 


tiply the expression for ¥(r) as given by Eq. (7) _ 


by the corresponding expression for ¥*(r) with 
r set equal to fo, where the asterisk denotes the 
complex conjugate of the quantity to which it 


’ Whittaker and Watson, Modern Analysis (Cambridge 
University Press, 1927), Chap. XI. 
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is attached. This yields 
—Yo(r) ] 
+2 (re) * E*(to, ry) 


The prime on the last summation indicates that 
terms for which j’=j’’ are to be omitted from 
the sum. If we now take the configurational 
average of both sides of this equation we obtain 


(ro) 
—Yo*(To) ] 
—Yo(r) ] 


+ f H(ry) 
Vv 


X E(x, ry) E*(ro, ry)dry 
N-1 


x ry) E*(fo, ty )drydr;, (19) 
where 


H(t) = f (20) 


We now make a similar assumption to that made 
in the preceding section, although the physical 
meaning is not so clear in this case. We assume 
that it is valid to replace 


Making this substitution and replacing (V—1)/N 
by unity since we are assuming JN is large (com- 
pared to unity) we obtain the integral equation 


(To) 
—Yo*(ro) 
—volr)] 
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By substituting the right-hand side of Eq. (9) 
for ((r)) and the complex conjugate for 
(y*(fo)) in the above equation we obtain a 
simpler form of the equation: 


Vv 


Now the above integral equation for ) 
is in too complicated a form to obtain practical 
solutions of it directly. However, we can simplify 
the equation further in such a way that reason- 
ably rapidly convergent series solutions may be 
obtained in many cases. An examination of the 
integral equation (22) shows that by obtaining 
the Liouville-Neumann iteration solution of the 
equation, considered as an integral equation in 
— W(t) ] with the first term 
on the right in Eq. (22) considered as the in- 
homogeneous term of -the integral equation, 
(v(r)¥*(r0)) satisfies the integral equation 


= 


(22) 


where the kernel L(r, 1; 1’) is given by the 
iteration series: 


Lit, fo; = M,(f, fo; fr’), 


p=0 


(24) 


M,(t, fo; = f f 
X r’)E(r, r’’) 
XE*(ro, 
to; r') =H E(x, r’)E*(r0, r’). 


(25) 
(26) 


This solution has the disadvantage that the 
series solution for the kernel L(r, ro; r’) will be 
very slowly convergent since the factors E(r, r’’) 
and E*(ro,r’’’) in the iteration integrals are 
imaginary exponentials and do not promote con- 
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vergence of the series, in general. An improve- 
ment on the situation is possible by obtaininga . 
more rapidly convergent series for the kernel. 
To do this we first find the partial differential 
equation satisfied by L(r, ro; r’). lf we apply the 
operator (Vo?+o")(V?+ ko") to both sides of Eq. 
(24), where V? is the Laplacian operator in the r 
coordinates and V,? the Laplacian operator in the 
ro coordinates, the result of this operator acting 
on Mp, will be while the 
result of the operator acting on the remainder 
of the series is readily found to be 162°G(r) - 
XG*(ro)L(r, ro; r’). Thus the integral equation 
kernel satisfies the partial differential equation 


— 164°G(r)G*(r0) JL(r, ro; r’) 
(27) 
If we now add 

{4xG(r)[Vo? +k**(r0) ] 
to; 1’) (28) 
to both sides of this equation, and remember that 
k?(r) =ko?+4xG(r), (16) 
= ko? +44G*(t0), (16’) 


and 


the equation can be written 
= 
+ {4rG(r)[Vo? +k**(r0) 


+4nG*(ro) [V2 +k2(r) ]}L(r, ro; 7’). (29) 


Let us now turn our attention to the function 
K(r, r’) defined as the solution of the equation 


V*K(r, r’)+k2(r)K(r, —4n8(r—r’), (30) 


which represents outward travelling waves at 
infinity. It will be noted that K(r, r’) bears the 
same relationship to Eq. (15) as E(r, r’) does to 
Eq. (4). K(r, r’) is a Green’s function for Eq. 
(15) and can be interpreted as the wave field 
produced by a point source of unit strength 
placed at the point r’ in a medium in which the 
propagation constant at the point r is given by 
k(r). Then K(r, r’\K*(ro, r’) satisfies the equa- 
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tion 
[Vo?-+h**(ro) ]K(r, 2’) K*(ro, r’) 
(31) 


From this one readily finds that L/(r, ro; r’) 
satisfies the integral equation 


L(r, =K(r, r’)K*(ro, r’) 
1 


YL", 
XK(r, K* (ro, (32) 


as may be verified by substitution in Eq. (29). 
The Liouville-Neumann iteration procedure ap- 
plied to this integral equation gives us a new 
series solution for the kernel of our integral 
equation: 


p=0 


1 


XLpa(r”’, 2”; 2’) r’’) 
| XK*(ro, (34) 
Lo(r, to; = K(r, r’)K*(ro, r’). (35) 


This, in general, will be a quite rapidly converg- 
ing series since K(r,.r’) contains a real negative 
exponential. The successive series terms will 
contain increasingly larger numbers of integra- 
tions involving the negative exponential and will 
thus die out rapidly. Thus, for example, the 
second term in the series is 


Li(r, = — f G(r") K(e", 


XK(r, r")K*(fo, r’)dr”’ 


re f (36) 
Vv 
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The solution to our problem has thus been 
reduced to solving the integral equation 


(vr) |*)= | [2 
+ f *)L(x, (37) 
4 


In case the kernel L(r, r; r’) can be accurately 
represented by the first term in the series (33), 
the integral equation becomes 


|?)= | |? 


This latter equation is probably the most im- 
portant practical result of our investigation. 


AUXILIARY PHYSICAL QUANTITIES 


There are several other quantities of im- 
portance in wave propagation theory and we 
shall now obtain expressions for the configura- 
tional average of these. The first of these is the 
gradient of the wave function. It is readily ob- 
served that the operation of taking the gradient 
and the integrations involved in obtaining the 
configurational average commute so that the 
average of the gradient of the wave function is 
simply the gradient of the configurational aver- 
age of the wave function: 


(38) 


The other quantity of importance is the flux 
which is proportional to the quantity 


(40) 


We may obtain the configurational average of 
this quantity by taking the gradient of both 
sides of Eq. (23) with respect to r and subtracting 
from this the gradient of both sides of the same 
equation with respect to ro. On subsequently 
setting rfo=r, we obtain 


+ f H(e’)(| v(x’) 
*(r, ro; (41) 
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(41) 


Again if L(r, ro; r’) can be replaced by the first 
term in its iteration series, this result can be 
written 


= (2) — H(z) 


~K(r, r’)VK*(r, r’) (42) 


Thus (S) can be calculated when (|¥(r)|*) is 
known. 


PHYSICAL INTERPRETATION 


In order to interpret the results which we have 
obtained, it is useful first to consider the scat- 
tering of waves by a single scatterer. In this 
case the wave function becomes simply 


¥(r) =Yo(r) 


which is the sum of the incident wave and the 
spherical wave scattered from the scatterer. If 
we consider a plane incident wave 


Yo(r) =A exp [—iky-r], 


and integrate the flux per unit area in the scat- 
tered wave, where y, repre- 
sents the second term on the right in Eq. (43), 
over a sphere whose center is at the position of 
the scatterer, we readily find that this total flux is 


4riky | 21 | 242. (45) 


(43) 


(44) 


' If we divide by the flux per unit area normal to 


ky in the incident wave, ikoA*, we obtain what is 
usually referred to as the scattering cross section, 

4n | £1 | 3, 


(46) 


It represents an area which when multiplied by 
the incident flux per unit area, gives the total 
flux in the scattered wave. 

Now if one integrates the flux per unit area 
caused by both incident and scattered waves 
[that is y*Vy—yVy*, where y is given by (43) ] 
over the same sphere, one obtains 


(47) 


The negative of this quantity represents the net 
inward flux through the sphere. Dividing the 
latter by the flux per unit area in the incident 
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wave gives the so-called absorption or capture 
cross section, o., since it represents the area 
which when multiplied by the incident flux per 
unit area normal to Ko, gives the flux which is 
effectively lost from the wave field. Thus 


—4n|g1|?— 


(48) 
where 


—4nIm(g1)/ko (49) 


is usually referred to as the extinction cross 
section, being the sum of the capture and scat- 
tering cross sections. In practical problems, the 
imaginary part of g; is usually negative so that 
the extinction cross section is positive.’ Ex- 
pressed in terms of the more familiar cross 
sections, the scattering coefficient g: becomes 


Os Rove (50) 


We may then define a scattering cross section 
per unit volume, S,(r), by 


S.(1) = f 
f le(s)|*n(x, s)ds=4eH(r), (51) 


and an extinction cross section per unit volume, 
S.(r), by 


S.(r) = f o-(s)n(r, s)ds 


= 


(52) 


Let us now examine Eq. (15) which governs 
the propagation of the average value of the wave 
function in the case where G(r) is a constant 
independent of position. In this case the plane 
wave solution of the differential equation will be 


*If there is no “true” capture, o, must be zero. For a 
particular field, this can be'shown from the field equations 
and, where appropriate, the boundary conditions at the 
scatterer. However, £ will still be complex because of 
“radiation resistance. 
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of the form es 
(W(r))=A exp [—ik-r], (53) 


where |k| =& and the direction of k is the direc- 
tion of propagation of the wave. Now since k is 
in general complex the wave represented in Eq. 
(53) is attenuated with distance. Now we have 


(54) 


so that when the number of scatterers per unit 
volume is sufficiently small we may expand the 
square root for k to obtain - 


21 
k~~ky+—Re(G) +—Im(G). (55) 
ko ko 


Thus the flux in the wave (53) is reduced by the 
factor 


2a 3 
[ex az Se (56) 


per unit distance. Hence in this case the medium 
acts as if it removes in a unit distance a fraction 
equal to the extinction cross section per unit 
volume of the flux carried by the wave. Alter- 
natively, we may say that the medium acts as 
if each scatterer removes from the wave an 
amount of flux equal to its extinction cross 
section times the flux per unit area incident on 
it in the average wave. This flux is partially ab- 
sorbed (if the capture cross section is not zero) 
and partially scattered out of the average wave 
propagating through the medium. It should be 
noted that this simple picture is applicable only 
when the number of scatterers per unit volume 
is sufficiently small so that the above expansion 
of the square root is valid. Otherwise there is 
interference between the scatterers and they no 
longer can be considered to absorb and scatter 
independently. 

In addition the phase velocity c(r) of the 
average wave is different from the phase velocity 
in the scatterer-free medium: 


1 
(57) 


Co? w? 


This may be attributed to the fact that there is 
a difference in phase between the wave incident 
on a scatterer and the scattered wave from the 


scatterer. When the scattered waves are com. 
bined with the incident wave there is a pro. 
gressive advancement or retardation of the phase 
of the disturbance and the result appears as a 
change in the phase velocity.’ 

Let us now turn our attention to the physical 
interpretation of our integral equation for the 
configurational average of the square of the 
wave function written in terms of the scattering 
cross section per unit volume [see Eqs. (51), 


(37) J: 
|?) = | |? 


1 
+— f (58 
4n 


Assume first that we may replace L(r, r; r’) by 
the first term in the series (33), | K(r, r’)|*. Now 


we have already pointed out that K(r, r’) can’ 


be interpreted as the value of the wave function 
at the point r caused by a unit point source at 
the point r’, in the continuous medium [charac- 
terized by the propagation constant k(r)] in 
which (¥(r)) is propagated. Thus in this case 
the integral equation above simply states that the 
value of {|¥(r)|?) at any point r is equal to the 
mean-square value of the average wave at the 
point r, |(¥(r))|?, plus contributions owing to 
scattered waves from all other points of the 
medium, the strength of the sources per unit 
volume at these other points being proportional 
to the value of (|¥(r)|*) at these points and to 
the scattering cross section per unit volume at 
these points. The various contributions, it will 
be noticed, combine as if the waves to which they 
are owing are incoherent, that is, bear a random 
phase relationship to one another. If we obtain 
the series solution of this integral equation by the 
Liouville-Neumann method, this solution, too, 
has a simple physical interpretation. The solution 
can be written 


(v(x) = Vals), (59) 
m=0 
1 
Wo(r) = | (61) 
7 The condition for no dispersion in the of the 


scatterers is obviously: g(w, s)=w*y(s), where y is inde 


pendent_of w. 
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Thus the first term in the series represents the 
contribution caused by the average wave, the 
second term, the contributions caused by a single 
scattering of the average wave from each point 
in the medium to the point in question, the third 
term, the contributions caused by scattering 
twice of the average wave, etc. Hence, when 


L(t, r; t’) can be replaced by the first term in - 


its series solution, we have a very simple and 
logical physical interpretation of the integral 
equation which we have derived. 

When the full series solution for L(r, r; r’) is 
substituted in the integral equation, however, 
the meaning of the additional terms is not so 
clear. An examination of them indicates that 
they are some sort of residual interference terms 
which are not eliminated in the averaging process. 
This is apparent from the fact that the imaginary 
exponential terms are still present [Eqs. (34), 
(36) ]. The physical interpretation of each term 
is quite simple: Thus, the terms corresponding to 
the next term in the series for L, that is ZL; [Eq. 
(36)], may be interpreted as the resultant at 
the point r of a wave scattered from the point r’ 
to r by two paths, a direct path from r’ to r and 
also by propagation from r’ to another point r”’ 
from which the wave is rescattered coherently to 
the point r, the result being summed over all 
intermediate points r”’. Similarly the higher 
terms of the series correspond to larger numbers 
of intermediate coherent scatterings. It is dif- 
ficult to see, however, why these additional 
terms should be present in the equation. A rough 
calculation of the order of magnitude of these 
terms shows that in many practical cases, they 
are not very significant. Thus the ratio of the 
terms caused by L, to the terms caused by 
Ly=|K(r, r’) |? seems to be of the order of mag- 
nitude of (k?—ko*)/ko’?, so that if the change in 
the velocity of the waves in the medium pro- 


duced by the introduction of the scatterers is 


small, these higher terms are of little importance. 

.In conclusion, it may be helpful to point out 
that in most measurements on waves it is the 
mean-square value of the wave function rather 
than the wave function itself that is measured; 
therefore (|¥(r)|*) will usually be of more inter- 
est than (¥(r)). Now the difference between 
(|¥(r) |?) and |((r))|? can be interpreted from 
the view of our statistical ensemble of con- 
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figurations as effected by the statistical fluctu- 
ations of ¥(r) about (¥(r)) in different members 
of the ensemble. Another method of viewing this 
difference may be helpful in bringing out the 
physical interpretation. Consider instead of a 
statistical ensemble of configurations of scat- 
terers, a collection of scatterers which are moving 
about slowly (so that Doppler effects may be 
neglected) in such a way that over a period of 
time the collection passes through all the con- 
figurations of our original ensemble. In this case 
{|¥(r)|*) can be interpreted as the long time 
average of |y/(r, ¢)|*, being the value of 
the wave function at the point r at the time ¢. 
For. such a collection of scatterers (¥(r)) has 
little physical significance unless 


To define the meaning of (¥(r)) precisely, one 
may say that if y(r, 4) represents the value of 
the wave function at the point r at time ¢, then 
(¥(r)) is the function which makes the integral 


five (63) 


a minimum, the integration being taken over a 
length of time sufficient for the collection of scat- 
terers to have passed through all the con- 
figurations of the ensemble. 

On the other hand, we have still to indicate 
why in scattering by a fixed configuration of 
scatterers, ‘‘incoherent scattering” is observed. 
The answer to this question lies in the fact that 
the regularly propagated wave or “coherent scat- 
tering” is precisely defined only as an average 
over the ensemble of configurations. For a fixed 
configuration, all of the scattering is strictly 
coherent but is artificially divided into a part 
which is estimated to propagate uniformly and 
a remainder which is referred to as ‘incoherently 
scattered.”” The basis of this distinction lies in 
the fact that for a collection of a large number 
of scatterers a particular ‘‘unprepared” collection 
will have particular physical properties which 
do not deviate greatly from the average physical 
properties of a properly defined statistical en- 
semble of collections because of the lack of “‘cor- 
relation” implied in the word “unprepared”’ as 
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to the positions of the individual scatterers. This 
allows one to estimate the properties of the aver- 
age over an ensemble from an experiment on a 
particular configuration with a very high prob- 
ability that the estimate is correct, although the 
possibility exists of a wide deviation if the se- 
lected ensemble should be particularly strongly 
ordered. In experiments usually several samples 
are used so as to guard against this remote pos- 
sibility. We cannot enter further into this dis- 
cussion of the relationship between statistical 
concepts and experiment, but refer the reader to 
almost any text on statistical mechanics for 
further elucidation. 


DISCUSSION 


It is of interest to indicate the directions in 
which the present theory may be generalized to 
embrace a wider range of phenomena: 

(1) We may attempt to lift the restriction to 
isotropic (s-) scattering. This could be accom- 
plished by including higher order scattering 
terms. Thus, in the next approximation we would 
include terms corresponding to dipole (p-) scat- 
tering which are of the form 


B-VE(r, r;), (64) 
where B is some vector function of the external 
field acting on the scatterer, such as its gradient. 
By adding further terms corresponding to quad- 
ripole scattering, etc., one could, in principle, 
adequately treat any type of scattered wave. The 
difficulties in this extension of the theory appear 
to be simply ones of increasing complexity of the 
expressions handled. 

(2) One might attempt to consider more 
general forms of probability distribution func- 
tions than those included in Eq. (1). Such more 
general functions would involve correlations 
between the positions of various scatterers so 
that the probability of finding a particular scat- 
terer in the neighborhood of some point would 
be a function of the positions of one or more of 
the other scatterers. This type of function would 
allow a more accurate representation of the 
state of affairs in ‘‘amorphous’”’ solids, liquids, 
and imperfect gases where such correlations do 
exist. Tempting as such a generalization might 
be, the well-known difficulties which arise in 
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averaging over correlative probability distriby- 
tion functions, are discouraging of elegant 
attempts at a rigorous theory in this direction, 
Some promise might appear, however, in an 
approximate treatment of this problem by sim- 
plifying the correlation between scatterers to 
that between the position of each scatterer and 
the average distribution of the remainder such 
as is done in the well-known approximate treat. 
ments of cooperative phenomena in statistical 
mechanics, or by the use of such methods as are 
employed in the theory of x-ray (single) scatter. 
ing in liquids. 

(3) A further generalization may be made in 
the extension of the above treatment to vector 
waves. If the vector wave function satisfies the 
vector wave equation and if the wave scattered 
by a scatterer is proportional directly to the 


external field acting on the scatterer, then the 


generalization follows immediately. However, in 
the case of greatest practical importance, elec- 
tromagnetic waves, the above is not true. The 
wave scattered by a scatterer is to first order 
proportional, not to the external field acting on 
the scatterer, but to the curlcurl of this field! 
This introduces singularities into the integrands 
of the integral equations which must be omitted 
from the integration by enclosing them by small 
spheres. These singularities lead in the treatment 
of the average value of the wave function to the 
well-known Lorentz-Lorenz effect,* but the deri- 
vation along rigorous lines is long and involved. 
These singularities appear to cause still more 
difficulty in the treatment of the average value 
of the square of the wave function (electric or 
magnetic field strength) but there seems no 
reason why the analysis cannot be carried 
through. The problem is of sufficient practical 
importance as well as theoretical importance to 
warrant additional study. ‘4 
(5) Another generalization could be attempted 
in the direction of considering changes of fre- 
quency of scattered waves produced through the 
Doppler effect by moving scatterers. This ‘is 
closely related to the problem of inelastic scat- 


8 The vector wave function is the Hertz vector Z in 
this case; the scattered waves are of the form pE(r, %), 
where p, the dipole moment of the scatterer is pro : 
to the external electric field acting on the scatterer, which 
is in turn proportional to the curlcurl of the external Z 
field acting on the scatterer. 
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tering of light and fundamental particles by free 
or bound atoms in the quantum theory, since a 
change of frequency of the wave function cor- 
responds in quantum mechanics to a change in 
energy of the associated particle. Such a theory 
would be developed along lines in some ways 
similar to the treatment of diffuse spots in the 
theory of x-ray diffraction which are also the 
effect of inelastic scattering of photons by the 
crystal lattice. Such a theory might have con- 
‘siderable value in investigating experimentally 
the dynamics of thermal agitation in amorphous 
solids, liquids, and gases. 
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(6) Finally in a similar manner, the theory 
could be extended to include non-steady state 
scattering, in which the incident wave function 
does not have a harmonic dependence on the 
time. 
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Letters to the Editor 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 


Double Bragg Reflections of X-Rays 
in a Single Crystal 
C. J. Davisson anD F, E, HawortH 


Bell Telephone Laboratories, Murray Hill, New Jersey 
January 25, 1945 


T has been called to our attention that experiments 
similar to those which we recently described under the 
above title in these columns,' though made on a different 
substance, were made and published eight years ago by 
M. Renninger in Germany. 

It is a matter of some embarrassment that we were 
unacquainted with Renninger’s work, and of regret that we 
failed, in what was obviously an incomplete search of the 
literature, to find the articles in which it is described. 

The phenomenon which we observed in quartz and 
designated as “double Bragg reflection”’ is identical with 
that which Renninger?® observed in diamond and designated 
by the somewhat inappropriate name of Umweganregung, 
which may be translated as “circuitous excitation” or 
“indirect excitation.” 

1C, J. Davisson and F. E. Haworth, Phys. Rev. 


ane Zeits. f. Physik 106, 141 (1937); Zelte, f. Krist. 


The Red-Giant Stage of Stellar Evolution 


G. Gamow 
The George Washington University, Washington, D.C. 
January 24, 1945 


HEN a normally evolving main-sequence star 
(point-source model with the carbon-cycle energy 
production) consumes all hydrogen contained in its central 
regions, one must expect the reaction zone to move out- 
wards thus leading to the so-called “‘shell-source’’ model! 
consisting of the isothermal dehydrogenized core and the 
radiative envelope. 

It was previously suggested by the author’ that the 
progressive growth of the energy-producing shell must 
result in a large increase of stellar radius, thus bringing the 
star into the red-giant region of the Hertzsprung-Russell 
diagram. This point of view encountered, however, serious 
difficulties owing to the fact that the exact mathematical 
treatment of the isothermal core model* seemed to lead to 
the conclusion that no such model is possible if more than 


ten percent of total stellar mass is to be contained in the 
core. 

It must be indicated here that the above paradoxical] 
result is entirely owing to the arbitrary assumption that 
the isothermal core of the star must be governed by the 
ideal gas law, and disappears at once if we take into 
account the possibility of degeneracy near the center of 
the core. Indeed the impossibility of fitting in a massive 
isothermal core of ideal gas (polytrope m= ©) rests on 
the fact that above a certain density limit, the gas pressure 
in the center of the core is unable to support the weight of 
the core thus leading to an unlimited contraction. Physi- 
cally, however, such a contraction will be rapidly stopped 
because of the decreased compressibility of gas in the 


_ degenerated region. The stable state of a star with a core 


exceeding ten percent of the total mass, will thus consist 
of three different regions: (1) degenerated nucleus, (2) iso- 
thermal layer of ideal gas, and (3) radiative envelope. 
The size and the luminosity of the star will be completely 
determined by the condition of degeneracy at the boundary 
between (1) and (2), and the condition of energy-produc- 
tion (essentially 72-107 °K) at the boundary between 
(2) and (3). Preliminary calculations of such composite 
stellar model, based on fitting of partially degenerated 
isothermal solutions for the core to the centrally con- 
densed solution for the envelope, indicate that the growth 
of energy-producing shell may increase the radius of the 
star by a factor of few hundred, and lower its surface 
temperature by a factor two or three. 

Since more massive stars develop the shell at a much 
higher rate, the locus of shell-source stars in the Hertz- 
sprung-Russell diagram must be represented by a line 
branching off from the main-sequence in the direction of 
higher luminosities and lower surface temperatures. As- 
suming the age of the main bulk of stars to be about 
310° years, we find that this branching off must take 
place near +2 absolute magnitude. It seems very likely 
that the red-giant branch of the stellar population II 
(elliptic nebulae, nuclei of spiral nebulae, and globular 
clusters) recently discussed by Baade,‘ corresponds to the 
group of (pregalactic) stars passing through this particular 
stage of their evolution. This would lead us to a very 
interesting conclusion that the members of the stellar 
population I (spiral arms) are cosmologically considerably 
younger than those of the population II. 

It may be also expected that with the gradual approach 
of the energy-producing shell to the stellar surface, the 


_ star will begin to eject its outer layers into the surrounding 


space, thus passing rapidly into the so-called Wolf-Rayet 
stage of evolution. The resulting loss of a large part of 
stellar mass must lead finally to the unveiling of the white 
dwarf, which was gradually growing in the center of the 
star during the red-giant part of its evolution. This view 
may find some support in the fact that the known white 
dwarf stars, which are apparently deprived of any nuclear 
energy-source, possess internal temperatures of the order 
2X10? °K which corresponds to the carbon-cycle source. 
This could be indeed interpreted by saying that those 
dense, slowly-cooling bodies still retain the temperature of 
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the isothermal core of the red-giant stars in which they 
must have been formed. 

The results of more detailed calculations of the above 
“three-layer” stellar model, and also the discussién of 
possible astrophysical consequences will be published in 
due course. 

1G. Gamow, Phys. . 53, 595 (1938); Chas. Critchfield and G. 
Gamow, Astrophys. J. bo. 244 (1939). 

2G. Gamow, Phys. Rev. 65, 20 (1944). 

3S, Chandrasekhar and L. R. Henrich, pm ang | J. 94, 525 (1942); 


s, Chandrasekhar and M. Schoenberg, Astrophys. J. 96, 161 (1942). 
5 Ww. Baade, Astrophys. J. 100, 137 (1944). 


A Note on the Kepler Problem in a Space of 
Constant Negative Curvature 


L. INFELD AND A. SCHILD 
Department of Mathematics, University of Toronto, Toronto, Canada 
January 10, 1945 


CHRODINGER' recently solved the “Kepler problem” 
in a spherical or Einstein universe and obtained the 
interesting result that the energy spectrum is discrete 
everywhere. It is instructive to compare this spectrum 
with that in an “‘open” universe of constant negative curva- 
ture, which is in fact Milne’s universe.? As will be shown, 
the spectrum consists of a finite number of (mostly nega- 
tive) energy levels in addition to a continuous spectrum. 
To our knowledge, this is the first quantum mechanical 
problem to exhibit a finite number of discrete energy 
values. 
The line element of the hyperbolic 3-space may be 
written in the form: 


ds* = R°da*+ sinh? a(d6?+sin? 6d¢*). (1) 


Letting R->~, a->0 such that aR=r remains finite, this 
line element reduces to that of Euclidean space, in which 
r, 6, @ are the usual polar coordinates. Schrédinger’s equa- 
tion is taken in the form 


Ay + (2u/h?)(E— V)y =0, 


_ with the proper Laplacian belonging to (1). 
The potential energy V(a) of an electron in a central 
Coulomb field must satisfy Laplace's equation 


(2) 


We find 
V=—(Ze*/R)(coth a—1), (3) 
if we demand that as and V-——dZe*/r as 
R-> aR—r. 
Putting ¥y=o(a) in Eq. (2), where the Y; are 
normalized spherical harmonics, we aus the radial 
equation 


(¢/da) (sinh? ado /da)+(A+1—2v) sinh? a 
+2» sinh a cosh a ¢—l(i+1)o=0; 
where = 2uR?E/h?—1 and 
Introducing the ‘“‘density function" p=o sinh a, we can 
factorize Eq. (4): 
coth a—v/(l+1)+d/da} 
X {¢+1) coth a—v/(l+1) —d/da}p 
= 20+ (+1)? 0; 


(4) 


(S) 
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coth a—v/l—d/da} {l coth a—v/l+d/da}p 

(5’) 
The only possible discrete solutions are immediately ob- 
tained by the factorization method. We have, in the 
notation of Infeld’s: paper,! 


An = 
Pn" =C sinh"*? a exp (—va/(n+1)); 
= 4 
X(lcoth a—v/l+d/da)pn'; 1=0,1--+n, (7’) 
since sinh a~}e* for large a, Eq. (7) shows that if the 
normalization integral for p,” is to exist, we must have 


(n+1)?<v; n=0,1---mo (8) 
where np is the greatest integer satisfying the inequality. 
Thus we have a finite number of discrete energy levels: 
n=0,1---m. (9) 


It is interesting to note that the number of discrete 
levels no+1~v! =(R/a)!, where a is the radius of the first 
Bohr orbit of the hydrogen-like atom. Taking R to be of 
the order of 107° cm (the usual magnitude of the radius of 
the universe considered in relativistic cosmology), we find 
that no is a large number of order 10'*. The highest discrete 
energy level Eng lies between — 3h?/2uR? and h?/2uR* and 
may thus be positive or negative ; all others are negative. 

The factorization method breaks down for the continuous 
spectrum. Putting x =coth a, Eq. (4) transforms into 


n=0,1,2---, 


(6) 
(7) 


(10) 


(x?—1)? 


The range of x (for O¢€a<@) is from » to 1. The 
only solution of this Riemann P-equation, bounded for 
large x, is* 
o = — 1) 
F(1+/+4[—M— 
1+1+4[—+ 
2(l+1); 2/(@@+1)). (11) 


This solution is continuous everywhere in our interval 
(including ©), except possibly at x=1. For negative \ we 
regain the discrete spectrum discussed above. 

If \ is positive, —\¥=i\' is imaginary. Introducing a 
for x in (11) and replacing o by the “density function” p 
we have 


p=sinh'*! a exp (—(1+i'+J)a) 
X 
2(l+1); 2); s=2e-*sinha. (12) 
We demand that p be continuous and bounded for 
O0<a<,. We have seen that p is continuous, and it is 
clear that the factor multiplying the hypergeometric func- 
tion in (12) remains bounded for large a. It remains to be 
shown that the hypergeometric function is also bounded 
as a>, i.e., as 2-1. 
Putting z= 1, we find the ratio of successive terms of the 
hypergeometric series is of the form‘ 
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It follows® that the infinite hypergeometric series oscil- 
lates finitely. Abel’s theorem® then shows that the hyper- 


geometric function in (12) remains bounded as 21 or 


This establishes the existence of a continuous eigenvalue 
spectrum of all positive or 


E>h?/2uR*. (13) 


In order to obtain an indication of how a relativistic 
treatment might affect this problem, the authors repeated 
the calculations by use of the Gordon-Klein equation. 
The general character of the solution was found to be 
unchanged. 

The authors wish to thank Professor A. F. Stevenson for 
helpful discussion. 


1E. Schrédinger, Proc. Irish Acad. A46, 9 (1940). 
Rev. 59, 737 and A. F. Stevenson, Phys. 


ey 9, ay Eq. (9) (1938), and A. G, 


3 E. T. Whittaker and G. N. Watson, Modern Analysis (Cambridge, 
199m) Secs. 10.7, 10.71, 10.72 and Chaps. 14 
a, a. Bromwich, "An Introduction to the Theory of Infinite Series 
(The Macmillan Company, New York, 1926), Sec. 79, p. 241, example. 
5 Reference 4, Sec. 79, p. 241, rule (ii). 
® Reference 4, Sec. 51, p. 149. 


Classical Theory of the Point Electron 


Mario SCHONBERG 
Department of Physics, University of Sio Paulo, Sio Paulo, Brazil 
December 18, 1945 


HE field produced by a point charge can be split in 
two parts: 


(1) The attached field, Fat; 


(2) The radiated field, Feaa. 


We assume the total field to be the sum of the attached, 
the radiated and the external fields and that the sum of 
the attached and radiated fields gives the usual retarded 


field rat of the particle. 

The radiated field behaves as an external field and 
reacts on the emitting charge. The attached field does not 
act on the generating charge but acts on other charges. 
Therefore, each charge is under the action of the retarded 
fields of the others, but only under the radiated part of its 
own field. 

The particle of charge e and mass m has two kinds of 
momentum and energy: (1) the kinetic four vector 


= Gkins 


and (2) the four vector of the acceleration energy and 


ds? 


The four vector é arises from the interaction of the 
charge and its radiated field. Since the particle is under the 
action of the external field and its own radiated field, 
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the-equations of motion ought to be: 


If we apply the laws of conservation of energy and mo. 
mentum we get: 


d dx? dx” 
5; (Chin +Gue) — ats (2) 


because the rate of loss of energy and momentum is 


according to the well-known Larmor formulae. Comparing 
Eqs. (1) and (2) we see that: 


e dx” _ 
c args 


4 Px, 
ds ds ds" (3) 


Equation (3) is satisfied by taking: 


§(Feet— (4) 


From Eq. (4) it results that: 


=4(Foot-+ Feav)- (5) 


We define the stress-tensor T*” of the field in which there 
are m point charges in the following way: 


4nT,’ = + Foe) 


Equation (6) ensures the conservation of the energy and 
momentum of the system of the field and particles, and 
leads to a finite field energy which is not always positive, 
The energy and momentum of the field derived from the 
stress tensor 7*” are the components of a four vector, 
A detailed exposition of the theory will be given elsewhere, 

Our theory leads to the Lorentz-Dirac classical equations 
of motion and to a finite field energy without introducing 


any energies and momenta of non-electromagnetic nature ° 


besides the kinetic ones. There are no subtractions of 
infinite quantities. 


The Radiation Field of a Point Electron 


J. Lettre Lopes* MARIO SCHONBERG 
Department of Physics, University of Sio Paulo, Sdo Paulo, Brasil 
December 18, 1944 


HE radiation field of the electron has been defined in 
different ways. The most usual is to take it as the 
part of the retarded field of the charge that depends on the 
retarded acceleration and varies inversely with the retarded 
distance to the charge. Though this definition presents 
many obvious inconveniencies, it has the advantage of 
giving the right amount of the radiation emission. Dirac 
has proposed the following definition of the radiation field 
of a point charge: 


uy 
Frad, D= Fret— Faav, (1) 


elect 
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_ pF” being the antisymmetric tensor of the corresponding 


uy 
_—— field: Frad,p the tensor of the radiation 


field, Fee that of the retarded field and <.. the tensor of 
the advanced field. Dirac’s radiation field has many 
advantages: it is finite on the charge’s world line, it is 
relativistically invariant and it is a solution of the homo- 
geneous Maxwell equations. It has an inconvenience; it does 
not give the correct radiation losses. Indeed; it can be 
easily seen that: 


if we take into account Dirac’s expression of Frad, p on the 
charge’s world line: 


xt=x, x=y, x*=z, (4) 


ds* = gapdx*dx® = (dx®)*— (dx')? — (dx*)?—(dx*)?. (5) 


The first term in the right-hand side of Eq. (2) is twice the 
rate of change of the four vector of acceleration energy and 
momentum and does not give rise to losses, the second term 
corresponds to twice the loss of energy and momentum by 
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radiation given by the Larmor formulae. This shows that 
instead of Dirac’s Frag, p it is better to take as the radia: 


tion field of a point charge F-.a: 
Fesa=4(Feet— Faav)- (6) 


Dirac’s classical equations of motion’ for a point charge 
can be put in the form: 
dx 


exit Fe rad), (7) 


ds? 


Fext being the tensor of the external field, arising from other 
charges. Equation (7) shows that the charged particle is 
under the action of the total field and its own radiation 
field defined according to Eq. (6), which should be reason- 
ably expected since only the radiation field which gets 
separated from the particle can react on it. 

A consistent theory of the radiation field of the electron 
that will be given elsewhere by one of us shows that all the 
difficulties of the classical theory of the electron can be 
overcome without any subtraction processes of the kind 
involved in the theories of Dirac and Pryce.” 


1P. A. M. Dirac, Proc. Roy. Soc, 167, 148 (1938). 

2M. H. L. Pryce, Proc. Roy. Soc. 168, 389 (1938). 

* Now at Princeton University with a fellowship of the U. S. Depart- 
ment of State under the provisions of the aevensee for the Promotion 
of Inter-American Cultural Relations. This work was done in 1943 
when J. Leite Lape was at the University of S$’ Pa Paulo with a fellowship 
of the Zerrener Foundation. 
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